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ABSTRACT 


A  study  was  conducted  of  a  new  testing  technique  using  flow  shaping  that  together 
with  some  modifications  of  the  AEDC  16-ft  Transonic  Propulsion  Wind  Tunnel  will  provide 
the  capability  to  test  full-scale  inlet/engine  configurations  with  forebody  effects  at  high 
maneuvering  angles  at  transonic  velocities.  The  method  used  to  obtain  the  flow  simulation 
for  high  maneuvering  angles  utilized  auxiliary  flow  shaping  and  geometric  pitch.  An 
analytical  potential  flow  method  was  used  to  determine  the  configuration  of  devices 
necessary  to  produce  the  required  flow  fields;  these  devices  were  then  checked 
experimentally  to  verify  the  results.  The  experimental  results  compare  favorably  with  the 
predicted  flow  fields,  and  the  use  of  flow  shaping  devices  for  simulation  of  the  inlet  flow 
field  for  certain  inlet  configurations  is  promising  for  angles  of  attack  up  to  20  deg  with 
angles  of  yaw  up  to  5  deg. 
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Cp  Pressure  coefficient,  (pL  -  pj/q„ 

K#  Circumferential  distortion  factor 

Ml  Local  Mach  number 

AMl  Ml  -  M„, 

M„,  Free-stream  Mach  number 

N  Any  free-stream  Mach  number  0.9 

pL  Local  static  pressure,  psfa 

pt2  Compressor-face  local  total  pressure,  psfa 

p^  Free-stream  static  pressure,  psfa 

ptoo  Free-stream  total  pressure,  psfa 

pt2/pt_  Steady-state  inlet  total  pressure  recovery 
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Y  Horizontal  coordinate,  positive  as  indicated 

Z  Vertical  coordinate,  positive  down 

a  Angle  of  attack,  deg,  positive  up 

1 3  Angle  of  yaw,  deg,  sign  as  indicated 

e  Upwasli  angle,  deg,  positive  up 

A e  (eM„=0.9  '  eM„=N)  or  (eC  yl.+  Inlet/Enginc  '  ^Cyl.  Only),  deg 

d  Local  flow  vector,  (tan'1v/tan2e  +  tan2  a),  or  angle  around  cylinder,  deg 

(see  Fig.  30) 

Angle  of  top  and  bottom  test  section  walls,  deg 
a  Sidewash  angle,  deg,  positive  away  from  Cylinder  No.  1 

Aa  (0M„=O.9  ‘  0M„=N  )  °r  (°Cyl.+  tnlet/Enginc  *  °Cyl.  Only  )-  deg 

Local  flow  orientation,  deg,  tair1  (tan  o)/(tan  e)  (see  Fig.  30) 
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SECTION  I 
INTRODUCTION 

A  major  problem  area  in  the  development  programs  of  highly  maneuverable  aircraft 
is  the  integration  of  the  propulsion  unit  and  airframe  into  an  efficient  operational  system. 
Ground  test  facilities  provide  an  economical  means  of  determining  solutions  to  design 
problems,  but  the  degree  of  full-scale  flight  simulation  is  generally  much  less  than  desired. 
The  performance  requirements  of  many  of  these  aircraft  may  include  considerable 
maneuvering  at  angles  of  attack  up  to  25  deg  in  the  transonic  Mach  number  range  as 
shown  in  Fig.  la  (Appendix  I).  As  the  angle-of-attack  requirement  increases,  the  inlet 
recovery  decreases  and  the  engine  face  circumferential  distortion  factor  increases  as  shown 
in  Fig.  lb  (Ref.  I).  This  results  in  a  decrease  in  the  engine  performance  at  maneuvers 
where  high  performance  is  essential  for  an  air  superiority  aircraft. 

Existing  propulsion  wind  tunnels  are  not  capable  of  testing  the  full-scale  inlet/engine 
systems  at  angles  of  attack  above  10  to  12  deg  simply  because  of  physical  limits  of  the 
tunnels,  and,  at  present,  there  exists  no  capability  to  test  at  a  pitch  and  yaw  attitude 
simultaneously.  This  angle-of-attack  limitation  is  reduced  further  because  of  the  necessity 
of  simulating  the  aircraft  forebody.  The  degree  of  forebody  simulation  required  depends 
on  the  aircraft  configuration,  but  certain  wind  tunnel  tests  (Ref.  2)  have  shown  that  the 
fore  body  may  produce  differential  flow  angularity  across  the  inlet  as  high  as  15  deg.  This 
leaves  a  considerable  gap  between  desired  and  available  test  angle-of-attack  capability.  This 
gap  is  shown  clearly  in  Fig.  la.  Placing  this  test  limitation  on  Fig.  lb  shows  that  the 
problems  of  reduced  inlet  recovery  and  increased  circumferential  distortion  are  only  starting 
to  occur  at  the  upper  limit  of  the  testing  capability  leaving  the  real  problem  to  be  solved 
during  flight  test.  Solving  the  inlet/engine  compatibility  problems  during  flight  test  can 
be  very  costly  as  demonstrated  clearly  in  a  recent  high  performance  aircraft  development 
program. 

Since  the  full-scale  engine  must  ultimately  operate  with  the  full-scale  inlet  in  the 
presence  of  the  airframe  throughout  the  required  operating  envelope,  it  is  highly  desirable 
that  the  full-scale  performance  be  simulated  in  a  ground  test  facility  before  flight  test. 
However,  using  present  testing  techniques,  this  simulation  would  require  a  new,  very  large 
propulsion  wind  tunnel. 

In  order  to  utilize  existing  test  facilities,  a  study  was  made  to  determine  new  test 
techniques  and  tunnel  modification  requirements  that  would  increase  the  range  of  flight 
conditions  for  testing  full-scale  inlet/engine  systems  in  the  AEDC  16-ft  Transonic  Wind 
Tunnel  (PWT-16T). 


SECTION  II 
OBJECTIVE 

The  objective  of  the  study  was  to  develop  a  testing  technique  and  to  produce  criteria 
for  the  modifications  of  the  AFDC  PWT-16T  that  together  would  provide  a  capability 
to  test  full-scale  inlet/engine  configurations  for  highly  maneuverable  aircraft  at  transonic 
velocities  over  an  angle-of-attack  range  up  to  25  deg  and  at  angles  of  yaw  up  to  5  deg. 
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SECTION  III 
FLOW  SHAPING 


3.1  INTRODUCTION 

The  new  testing  techniques  investigated  involved  the  generation  of  the  inlet  flow 
field  corresponding  to  a  maneuvering  condition  (that  is,  pitch  and  yaw)  through  the  use 
of  auxiliary  flow  shaping  (or  flow  deflection)  devices  in  the  vicinity  of  the  inlet.  The 
approach  taken  during  the  investigation  was  to  approximate  the  inlet  flow  field  for  the 
higher  angles  of  attack  through  a  combination  of  flow  shaping  and  geometric  positioning 
of  the  inlet/engine  system.  For  moderate  angles  of  attack  and  yaw,  the  desired  yaw 
component  of  the  flow  field  could  be  produced  by  flow  shaping  devices  and  the  upwash 
(or  downwash)  could  be  produced  by  either  flow  shaping,  geometric  positioning  of  the 
model,  or  a  combination  of  both. 

The  steps  taken  in  developing  the  new  test  technique  for  testing  full-scale  inlet/engine 
systems  at  simulated  pitch  and  yaw  were:  (1)  determine  for  a  given  aircraft  configuration 
the  flow  field  at  the  inlet,  (2)  develop  flow  shaping  devices  which  are  capable  of  producing 
the  desired  flow,  and  (3)  verify  the  ability  of  the  device  by  conducting  flow  survey  tests 
using  subscale  models  of  the  devices  in  the  AEDC  1-Ft  Transonic  Wind  Tunnel  (PWT-1T). 
The  success  of  such  an  approach  obviously  depends  on  the  ability  to  correctly  determine 
the  flow  field  produced  by  the  aircraft  and  the  flow  shaping  devices.  To  correctly  determine 
the  flow  fields,  both  the  one  to  be  simulated  and  the  one  produced  by  the  simulation 
device,  requires  the  use  of  either  experimental  data  obtained  from  wind  tunnel  tests  or 
predictions  made  by  analytical  procedures.  No  doubt,  using  only  a  wind  tunnel  and  a 
trial  and  error  method,  shaping  devices  could  be  designed  to  give  the  desired  flow  field. 
However,  this  would  be  both  costly  and  time  consuming  because  a  larger  number  of  models 
probably  would  need  to  be  fabricated  and  tested.  Therefore,  an  analytical  method  of 
determining  the  shaping  device  requirements  was  needed.  In  addition,  analytical  predictions 
of  the  aircraft  flow  field  could  eliminate  the  need  for  conducting  a  wind  tunnel  test  to 
survey  the  aircraft  flow  field. 

3.2  COMPUTER  ANALYSIS 
3.2.1  Analytical  Method 

The  analytical  method  used  to  calculate  the  flow  field  is  a  three-dimensional  potential 
flow  solution  which  utilizes  a  vortex  lattice  to  describe  the  model.  A  compressibility 
correction  using  Goethert's  Rule  has  been  incorporated  into  the  solution  for  use  at  high 
subsonic  Mach  numbers.  By  using  the  mathematical  model,  the  flow  angularity  and  local 
Mach  number  are  computed  at  the  desired  location  for  the  study  model.  Plots  of  flow 
streamlines,  constant  flow  angularity,  and  local  velocity  vectors  are  all  available  if  desired. 
Most  of  the  computer  analyses  were  made  at  Mach  number  0.9  since  this  was  determined 
to  be  near  the  upper  limit  of  Mach  number  for  this  analytical  method.  A  complete 
description  of  the  analytical  method  including  the  theory,  a  program  listing,  and  a  sample 
problem  is  given  in  Appendix  II. 
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3.2.2  Study  Configuration 

It  was  not  the  intent  during  this  study  to  determine  flow  simulation  for  a  specific 
aircraft  configuration;  however,  a  specific  configuration  was  analyzed  in  order  to  verify 
the  analytical  prediction  method  and  to  have  an  example  of  the  type  of  How  simulation 
required. 

Because  of  the  difficulties  encountered  in  trying  to  model  an  inlet  analytically,  it 
was  assumed  in  this  study  that  a  given  inlet  would  affect  a  given  flow  field  in  the  same 
manner  and  degree  regardless  of  how  the  flow  field  was  created.  Whether  or  not  this 
assumption  is  true  can  only  be  answered  by  testing  the  complete  vehicle,  but  the  assumption 
is  reasonable  and  the  simplification  was  introduced  into  the  analytical  calculations  because 
of  the  greatly  reduced  complexity  of  the  analytical  calculations. 

The  configuration  shown  in  Fig.  2,  which  is  a  fuselage  forebody  of  a  typcial  high 
performance  aircraft,  was  selected  for  the  study  because  wind  tunnel  data  of  flow  angularity 
surveys  from  previous  tests  in  the  AEDC  PWT-16T  (Ref.  2)  were  readily  available  for 
verification  of  the  analytical  method.  This  configuration  was  also  a  good  choice  because 
the  inlets  are  side  mounted  (typical  of  today's  high  performance  aircraft),  and  flow  surveys 
were  made  far  aft  where  the  forebody  effect  should  be  larger.  The  mathematical  model 
(location  of  vortex  filaments  and  velocity  control  points)  used  for  the  study  configuration 
is  shown  in  Fig.  3. 

The  analytical  data  for  one  model  position  are  shown  in  Fig.  4.  These  data  are  for 
an  angle  of  attack  of  25  deg  with  zero  yaw  angle  and  are  shown  as  lines  of  constant 
angles  of  upwash  and  sidewash.  A  comparison  between  the  analytical  data  and  wind  tunnel 
data  taken  at  the  model  position  indicated  in  Fig.  3  is  also  shown  in  Fig.  4  to  verify 
the  capability  of  predicting  the  desired  flow  field  with  the  three-dimensional  potential 
flow  solution  used.  The  figure  shows  the  data  to  agree  well  within  the  range  needed  for 
flow  simulation. 

A  preliminary  analysis  of  the  upwash  at  the  inlet  station  due  to  wing  lift  was  also 
made  to  determine  if  the  wing  effect  needed  to  be  simulated.  The  intersection  between 
the  leading  edge  of  the  wing  and  the  fuselage  was  at  Station  390  with  the  inlet  survey 
location  at  Station  356  (see  Fig.  3).  The  analysis  showed  the  wing  to  have  little  effect 
on  the  flow  angularity  at  the  inlet  station  for  Mach  numbers  of  0.9  and  0.75.  The  maximum 
differential  flow  angularity  was  less  than  two  degrees  in  upwash  and  less  than  one-half 
degree  in  sidewash  for  a  Mach  number  of  0.75  and  an  angle  of  attack  of  25  deg.  one 
degree  in  upwash,  and  one-quarter  degree  in  sidewash  for  a  Mach  number  of  0.9  and 
an  angle  of  attack  of  25  deg. 

3.2.3  Dual  Hollow  Circular  Cylinders 

Two  main  criteria  for  the  shaping  devices  aside  from  producing  the  flow  field 
simulation  are  that  they  have  minimum  tunnel  blockage  and  generate  as  little  turbulence 
in  the  test  region  as  possible.  The  first  device  analyzed  that  showed  promise  was  a  pair 
of  hollow  circular  cylinders  placed  in  a  cross  flow.  During  the  theoretical  analysis  made 
on  this  device,  the  spacing  between  the  cylinders  and  the  angle  of  attack  of  each  cylinder 
was  varied  until  a  flow  field  was  produced  which  resembled  the  flow  field  produced  by 
the  study  model  configuration. 
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An  analysis  was  made  with  the  cylinders  in  a  tunnel  and  in  uniform  free  flow  to 
ensure  that  the  presence  of  the  tunnel  walls  would  not  greatly  modify  the  analytical  results 
at  least  up  to  Mach  numbers  of  0.9.  Figure  5  shows  the  mathematical  model  of  the  cylinders 
in  the  wind  tunnel.  Figures  6  and  7  show  the  upwash  and  sidewash  for  the  cylinders 
in  the  tunnel  and  in  uniform  free  flow.  This  comparison  shows  only  slight  changes  in 
the  data  between  the  two  conditions  analyzed. 

The  mathematical  model  of  the  dual  hollow  circular  cylinder  configuration  with  the 
cylinder  diameter  and  spacing  that  were  used  in  the  experimental  study  is  shown  in  Fig. 

8.  The  theoretical  values  of  constant  angles  of  upwash  and  sidewash  are  shown  in  Fig. 

9.  The  analysis  of  this  configuration  predicted  that  increases  in  the  angle  of  attack 
simulation  of  six  degrees  could  be  obtained. 

3.2.4  Modified  Hollow  Cylinders 

In  an  attempt  to  get  a  flow  field  with  an  upwash  that  was  more  nearly  uniform 
in  the  vertical  plane  rather  than  the  horizontal  plane,  as  was  the  case  of  the  dual  hollow 
circular  cylinders,  and  to  get  better  yaw  simulation,  a  pair  of  modified  hollow  cylinders 
was  analyzed. 

These  devices  were  flattened  on  the  side  toward  the  tunnel  centerline  by  the  width 
of  one  radius.  Figure  10  shows  a  computer  plot  of  the  mathematical  model  of  this 
configuration  (a  cross-sectional  view  of  the  device  is  shown  in  Fig.  11).  Theoretical 
calculations  were  made  on  this  configuration  with  variations  in  pitch  on  both  cylinders 
and  in  combination  of  pitch  and  yaw  on  one  cylinder.  The  increase  in  upwash  for  a 
given  set  of  pitch  angles  was  less  with  this  device  than  with  the  original  dual  hollow 
circular  cylinders;  however,  the  upwash  did  become  more  nearly  uniform  in  the  vertical 
plane  as  desired.  By  decreasing  the  spacing  between  the  cylinder,  an  overall  increase  was 
shown  which  resulted  in  an  increase  in  the  predicted  angle  of  attack  simulation  of  seven 
degrees. 

The  results  of  the  analysis  made  with  combination  of  pitch  and  yaw  on  one  of  the 
cylinders  showed  that  this  configuration  could  be  used  to  obtain  the  yaw  simulation 
requirements.  The  computer  plots  with  constant  angles  of  upwash  and  sidewash  for  a 
typical  pitch-yaw  combination  are  shown  in  Fig.  11. 

Although  the  primary  objective  of  this  study  was  to  test  shaping  devices  that  would 
give  a  degree  of  simulation  of  the  study  configuration,  it  was  felt  that  basic  information 
about  the  flow-shaping  characteristic  of  other  devices  w'as  needed  for  future  use  on  other 
types  of  aircraft  configurations.  As  a  result,  two  other  sets  of  flow  shaping  devices  were 
chosen  and  analyzed  as  possible  candidates.  The  two  sets  w'ere  pairs  of  turning  vanes  to 
be  placed  in  the  flow  ahead  of  the  inlet.  A  description  and  the  results  of  the  analysis 
of  the  two  sets  of  vanes  are  given  in  the  next  tw'o  sections  of  this  report. 

3.2.5  Large  Turning  Vanes 

To  obtain  general  information  for  possible  future  flow-shaping  requirements,  a  set 
of  large  turning  vanes  which  spanned  the  tunnel  were  analyzed.  The  end  effect  of  the 
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tunnel  walls  was  neglected,  and  the  vanes  were  assumed  to  be  twodimensional  for  the 
analysis.  To  get  uniform  flow  entering  the  vanes,  the  walls  in  the  plane  of  the  vanes 
were  also  assumed  to  be  twodimensional.  The  mathematical  model  of  this  pair  of  vanes 
is  shown  in  Fig.  12,  and  resulting  flow  angles  using  both  vanes  or  either  vane  are  shown 
in  Fig.  13.  Analytically  tills  device  appeared  to  give  very  uniform  flow  across  the 
trailing-edge  plane. 

3.2.6  Small  Turning  Vanes 

For  further  information,  a  set  of  small  turning  vanes  was  analyzed.  These  vanes  were 
designed  to  shape  the  tunnel  flow  locally.  The  mathematical  model  of  these  vanes  is  shown 
in  Fig.  14.  The  vanes  were  assumed  to  be  in  a  uniform  free  stream.  The  theoretical  values 
of  constant  angles  of  upwash  and  sidewash  are  shown  in  Fig.  15.  Data  for  only  half  of 
the  survey  plane  are  shown  since  the  theoretical  analysis  was  made  using  only  half  of 
the  vanes  and  the  symmetry  capability  of  the  computer  program.  These  data  show  an 
inflow  around  the  convex  vane  (No.  2)  and  an  outflow  around  the  concave  vane  (No. 
1). 

3.3  WIND  TUNNEL  EXPERIMENTS 

After  theoretically  determining  the  shape  of  some  devices  which  have  predicted  flow 
fields  that  can  be  used  to  simulate  an  inlet  flow  field,  subscale  model  wind  tunnel  tests 
were  conducted  to  verify  the  theoretical  analysis.  This  experimental  effort  was  carried 
out  in  the  AHDC  PWT-IT.  A  description  of  the  test  hardware  and  models  is  given  in 
the  following  sections. 

3.3.1  Wind  Tunnel  (AEDC  PWT-IT) 

The  AEDC  PWT-IT  is  a  continuous-flow  nonreturn,  transonic  wind  tunnel  equipped 
with  a  two-dimensional,  flexible  nozzle  and  a  plenum  evacuation  system.  The  test  section 
Mach  number  range  can  normally  be  varied  from  0.2  to  1.50.  Total  pressure  control  is 
not  available  in  the  AEDC-PWT-1T,  and  the  tunnel  is  operated  at  a  stilling  chamber  total 
pressure  of  about  2,850  psfa  with  a  ±  5-percent  variation  depending  on  tunnel  resistance 
and  ambient  conditions.  The  stagnation  temperature  can  be  varied  from  80  to  120°F  above 
ambient  temperature  as  necessary  to  prevent  moisture  condensation  in  the  test  region. 

The  general  arrangement  of  the  tunnel  and  its  associated  equipment  is  shown  in  Fig. 
16,  and  a  schematic  of  the  nozzle,  test  section,  and  wall  geometry  is  shown  in  Fig.  17. 
A  detailed  description  of  the  tunnel  and  its  capabilities  is  given  in  Ref.  3. 

3.3.2  Flow  Angularity  Probe 

A  flow  angularity  probe  was  used  to  survey  the  flow  field  at  the  engine  inlet  location. 
Details  of  the  probe  are  shown  in  Fig.  18.  This  probe  was  1/4  in.  in  diameter  and  was 
designed  to  survey  an  area  two  inches  on  either  side  of  the  center  of  the  tunnel  and 
approximately  two  inches  above  and  four  inches  below  the  center  of  the  tunnel.  The 
probe  was  calibrated  in  the  AEDC  PWT-IT  with  data  being  obtained  through  a  Mach 
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number  range  from  0.5  to  1.4  at  probe  angles  up  to  ±33  deg  in  both  pitch  and  yaw. 
The  flow  was  repeatable  to  ±0.1  deg,  and  the  measuring  accuracy  of  the  probe  was  within 
0.5  deg  in  flow  angularity  and  0.05  in  local  Mach  number. 

The  probe  mechanism  was  mounted  to  the  sector  roll  mechanism  (see  Fig.  17)  which 
allowed  the  probe  to  be  rolled  to  get  up  to  two  inches  offset  from  either  side  of  the 
tunnel  centerline  for  positioning  of  the  probe  tip  in  the  z-axis.  The  sector  was  then  set 
to  traverse  up  or  down  to  allow  positioning  of  the  probe  tip  in  the  y-axis.  Position  changes 
in  the  x-axis  could  only  be  made  manually  with  the  tunnel  shut  down. 

3.3.3  Instrumentation 

The  tunnel  plenum  chamber  pressure  was  measured  with  a  self-balancing  precision 
manometer.  All  of  the  other  pressures  were  measured  with  differential  transducers 
referenced  to  the  tunnel  plenum  pressure  and  displayed  on  electromanometers. 

The  tunnel  stagnation  temperature  was  measured  with  an  iron-constantan 
thermocouple  and  displayed  on  an  indicating  potentiometer-type  recorder. 

3.3.4  Experimental  Models 

3.3.4. 1  Dual  Hollow  Circular  Cylinders  (Dual  Cylinders) 

Based  on  the  unit  size,  unit  space  between  the  cylinders,  the  size  requirement  and 
location  of  a  typical  engine  inlet,  which  were  determined  from  the  analytical  method, 
a  pair  of  hollow  half-circular  cylinders  was  designed  for  testing  in  the  AEDC  PWT-1T. 
Because  of  the  relative  large  cylinder  radius  to  separation  distance  required,  it  was  necessary 
to  use  half  cylinders  to  allow  space  for  a  normal  inlet/engine  installation.  The  cylinders 
could  be  individually  positioned  to  vary  the  spacing  between  them,  and  both  could  be 
pitched  in  5-deg  increments  up  to  40  deg.  In  addition,  one  of  the  cylinders  could  be 
yawed  up  to  10  deg. 

Figure  19  shows  a  schematic  of  the  general  layout  of  the  cylinders  installed  in  the 
tunnel,  and  Fig.  20  shows  an  installation  photograph.  A  sketch  showing  the  position  of 
the  probe  tip  for  the  survey  grid  and  sign  convention  used  with  the  dual  hollow  cylinders 
is  presented  in  Fig.  21.  The  x-axis  location  of  the  probe  tip  was  at  the  middle  of  the 
cylinders  for  the  majority  of  the  data  surveys. 

Five  pressure  orifices  were  located  around  Cylinder  No.  I  at  30-deg  intervals  6  in. 
from  the  leading  edge  of  the  cylinder  (see  Fig.  19).  Three  static  orifices  and  a  total  pressure 
probe  were  located  inside  the  cylinder. 

3.3.4.2  Modified  Hollow  Cylinders  (Modified  Cylinders) 

Based  on  information  determined  from  the  analytical  method,  a  pair  of  modified 
hollow  cylinders  was  designed  for  testing  in  the  AEDC  PWT-1T.  These  cylinders  were 
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the  same  basic  size  as  the  dual  hollow-circular  cylinders  but  were  flattened  on  the  sides 
by  a  straight  segment  equal  in  length  to  one  cylinder  radius.  These  modified  cylinders 
were  designed  to  fit  on  the  same  strut  system  used  by  the  dual  cylinders  and  could  be 
positioned  over  the  same  pitch  and  yaw  range. 

Figure  22  shows  an  installation  photograph  of  the  modified  cylinders.  A  sketch 
showing  the  position  of  the  probe  tip  for  the  survey  grid  and  sign  convention  used  with 
the  modified  hollow  cylinders  is  presented  in  Fig.  23.  The  x-axis  location  of  the  probe 
tip  was  at  the  middle  of  the  cylinders  for  the  majority  of  the  data  surveys. 

Five  pressure  orifices  were  located  around  Cylinder  No.  1.  Three  orifices  were  located 
on  the  flattened  side  and  two  on  the  leeward  curve  at  30-deg  increments  6  in.  from  leading 
edge  (see  Fig.  23).  Three  static  orifices  and  a  total  pressure  probe  were  located  inside 
the  cylinder  (similar  to  the  arrangement  in  Fig.  19). 

3.3. 4.3  Large  Turning  Vanes 

Using  the  contours  and  spacing  determined  by  the  analytical  method,  a  set  of  large 
turning  vanes  was  designed  for  testing  in  the  AEDC  PWT-1T.  These  vanes  were  strut 
mounted  to  the  tunnel  frame  and  spanned  the  tunnel. 

Figure  24  shows  a  schematic  of  the  general  layout  of  the  large  turning  vane  installation, 
and  Fig.  25  shows  an  installation  photograph.  The  survey  and  sign  convention  used  with 
the  large  turning  vanes  are  shown  in  the  sketch  in  Fig.  26.  The  x-axis  location  of  the 
probe  tip  was  located  at  the  trailing  edge  of  vane  No.  2.  This  was  the  same  plane  for 
which  the  theoretical  data  were  computed. 

Four  pressure  orifices  were  located  along  the  centerline  of  the  top  vane  at  2-in. 
intervals  starting  8  in.  from  the  leading  edge.  Five  pressure  orifices  were  located  along 
the  centerline  of  the  bottom  vane  at  2-in.  intervals  starting  5.5  in.  from  the  leading  edge 
(see  Fig.  24). 

3.3. 4.4  Small  Turning  Vanes 

A  set  of  small  turning  vanes  was  designed  which  used  the  contours  determined  from 
the  analytical  program.  These  vanes  were  mounted  to  the  struts  used  by  the  dual  hollow 
cylinders. 

Figure  27  shows  a  schematic  of  the  location  of  the  vanes  in  the  tunnel,  and  Fig. 
28  shows  an  installation  photograph.  The  location  of  the  probe  tip  for  the  survey  grid 
and  the  sign  convention  used  with  the  small  turning  vanes  are  shown  in  Fig.  29.  The 
x-axis  location  for  both  the  theoretical  analysis  and  the  experimental  surveys  was  at  the 
trailing  edge  of  the  vanes. 

3.4  RESULTS  AND  DISCUSSION 

Since  it  is  difficult  to  present  flow  angularity  data  for  three-dimensional  flow  with 
a  two-dimensional  plot,  the  majority  of  the  data  is  presented  two  ways.  First,  a  local 
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flow  vector  is  used  where  the  length  of  the  vector  indicates  the  magnitude  of  the  angle 
0,  which  is  the  resultant  local  flow  angle  relative  to  the  free  stream  and  is  essentially 
proportional  to  the  cross  flow  velocity,  and  the  orientation  (i p)  indicates  the  component 
of  upwash  and  sidewash.  A  diagram  of  the  local  vector  is  shown  in  Fig.  30.  The  second 
way  shows  a  graphical  display  of  the  data  at  constant  angles  of  upwash  and  sidewash. 

3.4.1  Dual  Hollow  Circular  Cylinders 

During  the  wind  tunnel  test,  data  were  obtained  at  Mach  numbers  of  0.6,  0.9,  1.0, 
and  1.2.  Cylinder  No.  1  was  positioned  from  0  to  40  deg  in  pitch  and  was  varied  in 
yaw  from  0  to  -10  deg  (nose  toward  tunnel  centerline).  Cylinder  No.  2  was  positioned 
from  0  to  30  deg  in  pitch  and  remained  at  0  deg  yaw  for  all  test  conditions.  Data  were 
obtained  at  three  different  x-stations  in  the  tunnel  over  the  y-z  survey  grid  as  shown 
in  Fig.  21. 

The  local  flow  vectors  showing  the  comparison  between  the  predicted  and  the 
experimental  flow  field  for  the  cylinder  configuration  indicated  are  shown  in  Fig.  31  for 
a  free-stream  Mach  number  of  0.9.  These  same  data  are  shown  in  Fig.  32  at  constant 
angles  of  upwash  and  sidewash.  In  Fig.  32,  the  wind  tunnel  survey  data  are  presented 
as  the  solid  circles  showing  the  location  of  the  probe  with  the  number  above  the  circles 
showing  the  measured  upwash  angle  and  the  number  below  the  circles  showing  the  measured 
sidewash  angle.  A  maximum  deviation  of  one  degree  is  noted  between  the  predicted  and 
measured  angles  at  any  one  probe  location. 

Figure  33  shows  a  comparison  between  predicted  and  measured  local  Mach  number 
for  the  same  test  condition.  Here  the  agreement  is  not  as  good.  However,  the  overall 
trends  are  the  same. 

As  indicated  in  the  model  description,  pressure  orifices  were  located  around  Cylinder 
No.  1  to  determine  if  any  anomalies  are  present  when  the  two  cylinders  are  varied  in 
the  proximity  of  each  other  that  would  cause  early  separation.  Figure  34  shows  a  plot 
of  the  pressure  coefficients  as  determined  from  the  static  pressure  measurements  for  pitch 
angles  of  30  and  40  deg  when  the  other  cylinder  is  pitched  to  20  and  30  deg,  respectively. 

Figure  34  indicates  that  the  flow  did  not  separate  with  30-deg  pitch  up  to  an  angle 
of  150  deg  around  the  cylinder,  whereas  the  flow  separated  between  120  and  150  deg 
around  the  cylinder  at  a  pitch  angle  of  40  deg.  By  using  Reynolds  number  based  on 
the  cross  flow  component  of  velocity  for  each  pitch  angle,  these  data  appear  to  he 
consistent  with  other  data  taken  with  a  cylinder  in  cross  flow,  and  the  presence  of  the 
other  cylinder  does  not  appear  to  create  any  premature  flow  separation  on  the  cylinder. 

Static  and  total  pressure  pickups  were  also  located  inside  the  cylinder  to  determine 
the  internal  flow.  The  data  from  these  locations  showed  that  the  flow  through  the  cylinders 
was  slightly  greater  than  free  stream  at  Mach  numbers  of  0.6  and  0.9.  However,  at  Mach 
numbers  above  0.9,  the  flow  through  the  cylinders  remained  constant  at  the  Mach  number 
0.9  value.  This  indicates  that  the  critical  mass  flow  ratio,  that  is  mass  flow  going  through 
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the  cylinder  to  mass  flow  capture  area,  must  occur  near  Mach  number  0.9  for  this  cylinder 
design. 

Figure  35  shows  the  variation  in  local  Mach  number,  upwash,  and  sidewash  with 
free-stream  Mach  number  at  four  probe  locations  when  Cylinder  No.  1  is  positioned  at 
30-deg  pitch  and  0-deg  yaw,  and  Cylinder  No.  2  is  positioned  at  20-deg  pitch  and  Q-deg 
yaw.  Figure  35a  shows  that  even  though  there  is  a  variation  in  local  Mach  number  from 
point  to  point  of  approximately  Mach  number  0.2  for  the  same  free-stream  Mach  number, 
this  variation  is  uniform  over  the  free-stream  Mach  number  range.  Therefore,  it  is  expected 
that  if  a  particular  local  Mach  number  is  needed  at  a  certain  location,  then  the  free-stream 
Mach  number  could  be  set  to  give  this  local  Mach  number. 

Figures  35b  and  c  show  the  delta  values  of  upwash  and  sidewash  over  the  free-stream 
Mach  number  range  as  compared  with  the  data  at  a  free-stream  Mach  number  of  0.9. 
This  figure  shows  a  maximum  deviation  of  approximately  2  deg  in  upwash  and  1.5  deg 
in  sidewash  over  the  free-stream  Mach  number  range.  This  indicates  that  if  a  change  is 
necessary  in  free-stream  Mach  number  to  get  a  correct  local  Mach  number,  the  flow  field 
simulation  will  still  be  within  the  desired  range. 

3.4.2  Modified  Hollow  Cylinders 

During  the  wind  tunnel  experiments,  data  were  obtained  at  Mach  numbers  of  0.6, 
0.9,  1.0,  and  1.1.  Both  cylinders  were  positioned  from  0  to  30  deg  in  pitch  and  Cylinder 
No.  1  was  varied  in  yaw  from  0  to  10  deg  (nose  away  from  tunnel  centerline).  Data 
were  obtained  at  two  different  x-stations  in  the  tunnel  over  the  y-z  survey  grid  as  shown 
in  Fig.  23. 

The  local  flow  vectors  showing  the  comparison  between  the  predicted  and  the 
experimental  flow  fields  at  a  free-stream  Mach  number  of  0.9  when  Cylinder  No.  1  was 
set  at  30-deg  pitch  and  0-deg  yaw  and  Cylinder  No.  2  was  set  at  20-deg  pitch  and  0-deg 
yaw  are  shown  in  Fig.  36.  Figure  37  shows  the  data  contours  for  the  same  conditions 
at  constant  angles  of  upwash  and  sidewash.  The  wind  tunnel  survey  data  are  presented 
in  Fig.  37  as  the  solid  circles  with  the  number  above  the  circles  showing  the  measured 
upwash  and  the  number  below  the  circles  showing  the  measured  sidewash.  The  maximum 
deviation  between  the  predicted  and  measured  angles  is  approximately  1.5  deg.  It  is 
interesting  to  note  that  the  analytical  method  underestimated  the  upwash  and  sidewash 
angles  for  the  modified  hollow  cylinders  and  overestimated  them  for  the  dual  hollow 
cylinders. 

Figure  38  shows  the  comparison  between  predicted  local  Mach  number  and  measured 
local  Mach  number.  Here  again  the  agreement  is  not  as  good  as  the  upwash  and  sidewash, 
but  the  trends  are  the  same. 

The  pressure  coefficients  for  the  ray  of  pressure  orifices  around  the  cylinder  are  shown 
in  Fig.  39  for  the  cylinder  pitched  30  deg  with  0-deg  yaw.  These  data  show  that  the 
flow  did  not  separate  up  to  an  angle  of  150  deg  around  the  cylinder. 
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Static  and  total  pressure  pickups  were  also  located  inside  these  cylinders  to  determine 
the  internal  flow.  The  data  from  these  pickups  showed  that  the  flow  through  the  cylinders 
was  slightly  greater  than  Mach  number  0.6  at  a  free-stream  Mach  number  of  0.6  but  did 
not  increase  as  free-stream  Mach  number  increased.  This  indicates  that  the  critical  mass 
flow  ratio  must  occur  near  Mach  number  0.6  for  the  cylinder  design.  The  cylinder  wall 
thickness  was  considerably  greater  for  the  modified  cylinders  than  for  the  dual  circular 
cylinders  (strictly  for  simplicity  of  fabrication)  which  reduced  the  internal  flow  area  thus 
lowering  the  Mach  number  at  which  the  critical  mass-flow  ratio  occurred.  This  caused 
the  modified  cylinders  to  block  the  tunnel  at  free-stream  Mach  numbers  above  1.1  and 
at  some  cylinder  pitch  angles  at  Mach  numbers  above  1.0.  Careful  design  of  the  cylinder 
should  increase  the  maximum  free-stream  Mach  number  at  which  the  cylinder  shape  could 
be  used. 

Figure  40  shows  the  variation  of  local  Mach  number,  upwash,  and  sidewash  with 
free-stream  Mach  number  at  four  probe  locations  when  Cylinder  No.  1  was  set  at  30-deg 
pitch  and  0-deg  yaw,  and  Cylinder  No.  2  was  set  at  20-deg  pitch  and  0-deg  yaw.  Figure 
40a  shows  a  variation  in  local  Mach  number  from  point  to  point  for  the  same  free-stream 
Mach  number  of  less  than  0.1  Mach  number,  and  the  variation  is  uniform  over  the 
free-stream  Mach  number  range.  This  should  allow  for  adjusting  the  free-stream  Mach 
number  to  set  any  local  Mach  number  within  the  range  of  the  data. 

Figures  40b  and  c  show  the  delta  values  of  upwash  and  sidewash  over  the  free-stream 
Mach  number  range  as  compared  with  the  data  at  a  free-stream  Mach  number  of  0.9. 
This  figure  shows  that,  at  free-stream  Mach  numbers  of  1.0  and  below,  the  maximum 
deviation  is  1.5  deg  in  upwash  and  0.75  deg  in  sidewash.  (From  Fig.  40a,  it  is  noted 
that  a  free-stream  Mach  number  of  1.0  gives  a  maximum  local  Mach  number  of  1.15.) 
At  a  free-stream  Mach  number  of  1.1,  the  data  show  one  point  to  deviate  up  to  1.75 
deg  in  upwash  and  2.5  deg  in  sidewash. 

Figure  41  is  the  flow  vectors  showing  the  comparison  between  the  predicted  and 
the  experimental  flow  field  at  a  free-stream  Mach  number  of  0.9  with  both  modified 
cylinders  set  at  30-deg  pitch  and  0-deg  yaw.  The  same  data  are  shown  in  Fig.  42  for 
constant  angles  of  upwash  and  sidewash.  The  wind  tunnel  data  are  again  shown  as  solid 
circles  with  the  number  above  the  circle  giving  the  measured  upwash  and  the  number 
below  the  circle  giving  the  measured  sidewash.  The  maximum  deviation  between  predicted 
and  measured  values  is  somewhat  higher  for  this  condition  resulting  in  a  maximum  upwash 
deviation  of  2.5  deg  and  a  maximum  sidewash  deviation  of  2  deg.  This  condition  gave 
the  maximum  area  over  which  there  was  a  uniform  upwash.  The  approximate  upwash  angle 
over  this  area  was  9  deg. 

Figure  43  shows  a  comparison  between  predicted  and  measured  local  Mach  number. 
Here  again,  the  trends  are  similar  as  before,  but  the  absolute  values  are  still  not  in  good 
agreement. 
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Figure  44  shows  the  variation  of  local  Mach  number,  upwash,  and  sidewash  with 
free-stream  Mach  number  at  four  probe  locations  when  both  cylinders  were  set  at  30-deg 
pitch  and  O-deg  yaw.  Data  were  only  obtained  up  to  a  free-stream  Mach  number  of  1 .0 
because  of  tunnel  blockage  at  these  settings.  The  deviation  in  local  Mach  number  with 
free-stream  Mach  number  is  shown  in  Fig.  44a  and  is  in  the  order  of  0.2  Mach  number. 

Figures  44b  and  c  show  the  delta  values  of  upwash  and  sidewash  over  the  free-stream 
Mach  number  range  as  compared  with  the  data  at  a  free-stream  Mach  number  of  0.9. 
The  deviation  in  upwash  over  the  free-stream  Mach  number  range  is  shown  to  be  0.2 
deg,  and  the  deviation  in  sidewash  over  the  fre^stream  Mach  number  range  is  shown  to 
be  0.75  deg. 

Figure  45  shows  the  effect  of  yawing  Cylinder  No.  1  at  a  constant  pitch  angle  of 
30  deg  with  Cylinder  No.  2  set  at  30-deg  pitch  and  0-deg  yaw.  Cylinder  No.  I  was 
positioned  nose  away  from  the  tunnel  centerline  at  0,  5,  and  10  deg.  Cylinder  No.  1 
was  yawed  to  create  the  flow  fields  necessary  for  inlet/engine  yaw  simulation.  The  data 
show  that  as  the  yaw  angle  is  increased,  the  sidewash  angle  increases.  It  should  be  noted 
that  when  the  cylinder  was  yawed  to  10  deg,  the  survey  area  was  decreased  because  of 
interference  between  the  survey  probe  and  the  aft  end  of  the  cylinder. 

3.4.3  Large  Turning  Vanes 

During  the  wind  tunnel  experiments,  data  were  obtained  at  Mach  numbers  of  0.6. 
0.9,  and  1.0  with  both  vanes  installed  and  at  Mach  numbers  of  0.6,  0.9,  1.0,  and  1.1 
with  vanes  individually  installed. 

The  local  flow  vectors  showing  the  comparison  between  the  predicted  and  measured 
flow  fields  at  a  free-stream  Mach  number  of  0.9  are  shown  in  Fig.  46.  These  same  data 
are  shown  in  Fig.  47  for  constant  sidewash  (or  upwash)  angles.  It  is  apparent  from  this 
figure  that  the  analytical  method  does  not  work  with  this  type  of  model.  It  should  be 
noted  that  although  the  vanes  span  the  tunnel,  the  flow  field  is  not  uniform  spanwise 
as  it  would  be  for  the  truly  two-dimensional  case  analyzed. 

Analysis  of  the  surface  pressure  data  shown  in  Fig.  48  indicates  that  the  flow  choked 
between  the  vane  at  Mach  numbers  above  1.0.  Data  for  the  individual  vane  tests  show 
the  vanes  to  act  similar  to  a  wing  surface  with  the  bottom  vane  having  a  shock  wave 
formed  on  the  upper  surface  at  Mach  numbers  above  1.0. 

3.4.4  Small  Turning  Vanes 

During  the  wind  tunnel  experiments,  data  were  obtained  at  Mach  numbers  of  0.6, 
0.9.  1.0,  and  1.1  with  both  vanes  installed  and  with  the  vanes  individually  installed. 

Figure  49  presents  the  local  flow  vectors  showing  the  comparison  between  the 
predicted  and  the  experimental  flow  fields  for  a  free-stream  Mach  number  of  0.9  with 
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the  small  vanes  located  four  inches  apart.  These  same  data  are  shown  in  Fig.  50  for  constant 
angles  of  upwash  and  sidewash.  In  Fig.  50,  the  wind  tunnel  survey  data  are  again  shown 
as  solid  circles  with  the  number  above  the  circle  indicating  the  upwash  angle  and  the 
number  below  the  circle  indicating  the  sidewash  angle. 

It  is  interesting  to  note  in  this  comparison  that  the  theoretical  method  predicted 
the  flow  trends  and  the  overall  values  quite  well.  The  deviation  between  the  theoretical 
and  measured  data  varies,  from  the  vane  that  is  convex  to  the  plane  of  the  How  toward 
the  vane  that  is  concave  to  the  plane  of  flow,  from  a  maximum  of  3  deg  on  one  side 
to  approximtely  5  deg  on  the  other  side.  Apparently  the  shorter  length,  outflow,  and 
the  lower  pressure  at  the  trailing  edge  created  by  the  design  of  the  convex  vane  tended 
to  eliminate  the  viscous  problems  that  were  present  with  the  large  turning  vanes.  Although 
the  flow  field  produced  by  the  small  vanes  does  not  appear  to  represent  the  required 
flow  field,  the  analysis  does  indicate  that  within  limits  the  theoretial  method  can  be  used 
to  predict  the  flow  field  for  this  type  of  model. 

3.4.5  Study  Configuration  Flow  Simulation 

The  final  comparisons  are  those  between  the  data  obtained  with  the  study  aircraft 
configuration  in  the  AEDC  PWT-16T  and  the  data  obtained  with  the  shaping  devices  in 
the  AEDC  PWT-1T. 

A  comparison  between  the  local  flow  vectors  obtained  at  Mach  number  0.9  with 
the  study  configuration  in  the  AEDC  PWT-16T  at  an  angle  of  attack  of  10  deg  and  a 
yaw  angle  of  zero  and  data  obtained  at  Mach  number  0.9  with  the  dual  hollow  circular 
cylinders  in  the  AEDC  PWT-1T  is  shown  in  Fig.  51.  Cylinder  No.  1  was  set  at  40-deg 
pitch  and  0-deg  yaw,  and  Cylinder  No.  2  was  set  at  30-deg  pitch  and  0-deg  yaw.  To 
get  the  simulation  shown,  it  was  necessary  to  locate  Cylinder  No.  1  as  indicated  in  the 
figure  and  to  adjust  the  data  for  a  geometric  pitch  of  the  inlet/engine  to  5  deg.  (To 
adjust  the  data  for  geometric  pitch,  the  assumed  geometric  pitch  angle  of  the  inlet/engine 
is  added  to  the  measured  upwash  angle.  The  new  upwash  angle  is  then  used  in  the 
computation  of  the  flow  vector.)  Thus,  the  forebody  effect  plus  a  net  gain  of  5  deg 
in  pitch  is  obtained  from  the  cylinders.  As  shown  in  Fig.  40,  if  the  local  Mach  number 
is  not  simulated,  the  free-stream  Mach  number  could  be  adjusted  to  give  the  correct  local 
Mach  number  without  greatly  affecting  the  upwash  and  sidewash  simulation. 

Figure  52  presents  the  local  flow  vectors  showing  the  comparison  between  data 
obtained  at  Mach  number  0.9  with  the  study  configuration  at  an  angle  of  attack  of  19.8 
deg  and  a  yaw  angle  of  zero  and  data  obtained  at  Mach  number  0.9  with  the  dual  hollow 
cylinders.  The  setting  of  the  dual  hollow  cylinders  is  the  same  as  in  Fig.  51.  However, 
to  get  correct  simulation  as  shown  in  Fig.  52,  it  was  necessary  to  move  Cylinder  No. 
1  nearer  to  the  inlet  and  to  adjust  the  data  for  a  geometric  pitch  of  the  inlet/engine 
of  12  deg.  This  is  the  maximum  geometric  pitch  available  with  the  full-scale  inlet/engine. 
Here  the  forebody  effect  plus  a  net  gain  of  7.8  deg  was  realized  with  the  flow  shaping. 
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Figures  53  and  54  present  the  comparison  between  the  experimental  data  obtained 
at  Mach  number  0.9  with  the  study  configuration  and  the  theoretical  results  computed 
at  Mach  number  0.9  for  the  modified  hollow  cylinders.  Figure  53  shows  the  flow  vector 
plot  for  the  study  configuration  at  an  angle  of  attack  of  10  deg  and  a  yaw  angle  of 
4  deg  (lee  side).  The  modified  hollow  cylinders  were  set  with  Cylinder  No.  1  at  30-deg 
pitch  and  10-deg  yaw  (nose  away  from  tunnel  centerline),  and  Cylinder  No.  2  at  20-deg 
pitch  and  0-deg  yaw.  The  theoretical  results  are  adjusted  for  a  geometric  pitch  of  the 
inlet/engine  of  10  deg.  Thus,  the  forebody  effect  and  yaw  angle  simulation  were  obtained 
using  the  flow  shaping  and  the  pitch  simulation  by  geometric  pitch. 

Figure  54  shows  the  same  data  compared  at  constant  angles  of  upwash  and  sidewash. 
From  this  figure,  the  actual  differences  in  the  numerical  value  of  the  angles  can  be  seen. 
It  is  noted  that  a  maximum  of  2  deg  difference  is  shown  between  the  experimental  data 
with  the  study  configuration  and  the  theoretical  results  with  the  modified  cylinders. 

Figure  55  shows  the  comparison  between  data  obtained  at  Mach  number  0.9  in  the 
AEDC  PWT-16T  with  the  study  configuration  and  data  obtained  at  Mach  number  0.9 
in  the  AEDC  PWT-1T  with  the  dual  modified  cylinders.  Here  again  the  study  configuration 
was  at  an  angle  of  attack  of  10  deg  with  an  angle  of  yaw  of  4  deg  (lee  side).  However, 
the  modified  hollow  cylinders  are  set  with  Cylinder  No.  1  at  30-deg  pitch  and  10-deg 
yaw  (nose  away  from  tunnel  centerline)  and  Cylinder  No.  2  at  20-deg  pitch.  The  difference 
between  these  data  and  the  previous  data  is  the  spacing  between  the  cylinders.  The 
theoretical  data  for  the  cylinder  are  for  a  spacing  of  1  54  in.  full  scale,  and  the  experimental 
cylinder  data  are  for  a  spacing  of  138  in.  full  scale.  In  Fig.  55,  the  modified  cylinder 
data  are  only  adjusted  for  a  geometric  pitch  of  the  inlet/engine  of  5  deg.  Thus,  the  forebody 
effect,  yaw  angle  simulation,  and  a  net  gain  in  pitch  angle  of  5  deg  are  obtained  from 
the  flow  shaping. 


SECTION  IV 
BLOCKAGE  STUDY 


4.1  INTRODUCTION 

A  simulation  technique  developed  using  shaping  devices  is  of  little  value  unless  there 
remains  sufficient  space  in  the  tunnel  for  the  inlet/engine,  and  the  tunnel  will  operate 
with  the  inlet/engine  installed.  Therefore,  throughout  the  design  of  the  shaping  devices, 
space  was  always  allowed  for  the  installation  of  an  inlet/engine.  The  final  phase  of  the 
study  was  to  determine  the  operating  characteristics  of  the  wind  tunnel  when  both  the 
shaping  devices  and  equivalent  full-scale  inlet/engine  are  installed  in  the  test  section.  Although 
the  majority  of  effort  in  this  phase  was  devoted  to  determining  wind  tunnel  operating 
parameters,  the  flow  angularity  probe  was  mounted  at  a  fixed  location  just  above  and 
in  front  of  the  inlet  to  determine,  if  possible,  the  effect  of  the  inlet  on  the  flow  angularity. 
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4.2  WIND  TUNNEL  EXPERIMENTS 

4.2.1  Wind  Tunnel  (AEDC  PWT-1T) 

The  AEDC  PWT-1T  is  described  in  Section  3.3.7  and  Ref.  3.  However,  attention 
should  be  called  to  the  plenum  evacuation  duct  and  scavenging  scoop  duct  shown  in  Fig. 
16  which  are  metered  to  determine  the  auxiliary  weight  flow  requirement  for  operation 
of  the  wind  tunnel.  Pressure  orifices  were  located  along  both  test  section  walls  and  in 
the  tunnel  diffuser  to  determine  the  Mach  number  variations  in  the  test  section  and  the 
tunnel  pressure  ratio,  respectively.  During  the  blockage  study,  the  top  and  bottom  test 
section  walls  were  diverged  one  deg,  and  the  diffuser  Haps  were  closed. 

4.2.2  Instrumentation 

The  basic  instrumentation  was  the  same  as  that  described  in  Section  3.3.3.  In  addition, 
however,  the  scavenging  scoop  and  plenum  airflows  were  measured  with  square-edged 
orifices.  The  general  location  of  the  metering  orifices  is  shown  in  Fig.  16.  The  pressure 
upstream  of  the  orifices  and  the  pressure  differentials  across  the  orifice  were  measured 
with  differential  transducers  and  were  displayed  on  electromanometers. 

The  orifice  airflow  temperatures  were  measured  with  iron  constantan  thermocouples 
and  were  displayed  on  indicating  potentiometer-type  recorders. 

4.2.3  Models 

The  shaping  device  models  used  in  the  blockage  study  were  the  dual  hollow  cylinders 
described  in  Section  3.3.4. 1  and  the  modified  hollow  cylinders  described  in  Section  3.3.4. 2. 

The  inlet/engine  model  used  was  a  one-sixteenth  scale,  two-dimensional,  supersonic 
inlet/engine  which  was  available  from  a  previous  wind  tunnel  blockage  study.  The  forebody 
was  removed  since  the  technique  is  to  simulate  forebody  effects  and  maneuvering 
conditions.  The  inlet/engine  model  could  be  manually  set  in  pitch  from  0  to  8  deg.  The 
flow  through  the  inlet/engine  could  be  varied  by  a  control  valve  in  the  scavenging  scoop 
line  to  simulate  different  engine  power  settings. 

Figure  56a  shows  an  installation  photograph  of  the  dual  hollow  cylinders  and  the 
inlet/engine  model  as  viewed  from  the  top  of  the  inlet.  Figure  56b  shows  the  dual  hollow 
cylinders  and  the  inlet/engine  model  viewed  from  the  wind  tunnel  nozzle  section.  In  both 
of  these  photographs,  Cylinder  No.  1  is  pitched  to  30  deg  with  0-deg  yaw  and  Cylinder 
No.  2  is  pitched  to  20  deg  with  0-deg  yaw. 

Figure  57a  shows  an  installation  photograph  of  the  modified  hollow  cylinders  and 
the  inlet/engine  model  as  viewed  from  the  top  of  the  inlet.  Here  Cylinder  No.  1  is  pitched 
30  deg  with  1 0-deg  yaw  (nose  away  from  tunnel  centerline),  and  Cylinder  No.  2  is  pitched 
20  deg  with  0-deg  yaw.  Figure  57b  shows  an  installation  photograph  of  the  modified 
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hollow  cylinders  and  the  inlet/engine  model  as  viewed  from  the  wind  tunnel  nozzle  section. 
Here  both  cylinders  are  pitched  30  deg  with  0-deg  yaw. 

4.3  RESULTS  AND  DISCUSSION 

4.3.1  Effect  of  Inlet/Engine  Installation  on  the  Flow  Angularity 

To  determine  how  the  installation  of  the  inlet/engine  would  affect  the  measured  flow 
field  using  only  the  dual  hollow  cylinders,  the  flow  angularity  probe  was  installed  with 
the  inlet/engine  model  during  the  blockage  study  testing.  The  probe  was  located  at  an 
x-station  and  y-z  axis  location  for  which  data  were  obtained  with  the  dual  hollow  cylinder 
flow  shaping  devices. 

Figure  58  shows  the  changes  in  the  upwash  and  sidewash  as  a  function  of  Mach 
number  that  resulted  from  two  settings  of  the  dual  hollow  cylinders  with  the  inlet/engine 
model  at  0-deg  pitch.  In  both  cases,  the  presence  of  the  inlet/engine  increased  both  the 
upwash  and  the  sidewash  angles.  With  Cylinder  No.  1  set  at  30-deg  pitch  and  0-deg  yaw 
and  Cylinder  No.  2  set  at  20-deg  pitch  and  0-deg  yaw,  the  upwash  was  increased 
approximately  1.5  deg  maximum,  and  the  sidewash  was  increased  approximately  2.25  deg 
maximum.  With  Cylinder  No.  1  set  at  40-deg  pitch  and  0-deg  yaw  and  Cylinder  No.  2 
set  at  30-deg  pitch  and  0-deg  yaw,  the  upwash  was  increased  approximately  2  deg 
maximum,  and  the  sidewash  was  increased  approximately  2.5  deg  maximum.  These  data 
have  the  expected  trends  in  both  cases  showing  the  maximum  increase  of  the  lower  Mach 
number  with  the  effect  decreasing  as  Mach  number  increases. 

4.3.2  Tunnel  Operating  Characteristics 

During  the  blockage  test,  tunnel  runs  were  made  for  most  positions  of  the  dual  hollow 
cylinders  and  modified  hollow  cylinders  previously  tested  with  inlet/engine  model  set  at 
0-  and  8-deg  angle  of  attack  with  simulated  engine  power  settings  of  idle  and  cruise.  It 
was  found,  for  the  most  part,  that  the  addition  of  the  inlet/engine  model  with  either 
pair  of  cylinders  had  little  effect  on  the  maximum  Mach  number  operating  capability 
of  the  tunnel  for  the  same  cylinder  settings  without  the  inlet/engine  model. 

During  the  blockage  test,  scavenging  and  auxiliary  weight  flows  required  for  operating 
the  tunnel  were  determined  by  metering  the  flow  through  the  scavenging  scoop  and  plenum 
evacuation  ducts.  Calculations  of  the  weight  flows,  measured  by  metering  orifices,  were 
consistent  with  the  procedures  outlined  in  Ref.  4. 

Figure  59  shows  the  ratio  of  the  auxiliary  weight  flow  to  the  tunnel  weight  flow 
required  to  operate  the  AEDC  PWT-1T  when  the  dual  hollow  cylinders  or  modified  hollow 
cylinders  are  installed  with  the  inlet/engine  model,  assuming  the  tunnel  is  started  with 
the  cylinders  and  inlet/engine  set  at  zero  pitch  and  yaw  until  after  the  engine  is  started. 
The  present  performance  curves  for  the  Plenum  Evacuation  System  (PWT-PES)  for 
operation  of  the  AEDC  PWT-16T  are  shown  in  the  figure. 
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By  assuming  that  the  AEDC  Engine  Test  Facility  airside  and  exhaust  plants  are  used 
lor  make-up  air  and  scavenging,  respectively,  thus  leaving  the  total  capacity  of  the  PWT-PES 
for  plenum  pumping,  it  appears  that  the  majority  of  the  requirement  for  operation  with 
the  dual  hollow  cylinders  and  inlet/engine  could  be  met.  However,  the  majority  of  the 
requirement  for  operation  with  the  modified  hollow  cylinders  and  the  inlet/engine  is  outside 
the  present  PWT-PES  capability. 

Two  additional  curves  are  shown  in  the  figure  which  represent  the  capacity  that  would 
be  available  with  the  addition  of  a  double  compressor  unit  equal  to  the  present  double 
units  or  the  addition  of  a  complete  third  increment  (five  units).  It  can  be  seen  that  the 
addition  of  a  double  unit  would  satisfy  the  requirements  for  operation  with  either  pair 
of  cylinders  with  the  inlet/engine  and  allow  some  leeway  for  How  shaping  devices  on 
inlet/engines  with  greater  tunnel  blockage. 

SECTION  V 

SUMMARY  AND  RECOMMENDATIONS 

The  analytical  and  experimental  data  obtained  during  the  How  shaping  study  with 
the  dual  hollow  circular  cylinders  and  modified  hollow  cylinders  show  excellent  agreement. 
The  analytical  and  experimental  data  obtained  with  the  small  turning  vanes  gave  fair 
agreement,  but  the  agreement  between  analytical  and  experimental  data  with  the  large 
turning  vanes  was  poor. 

Of  the  devices  tested,  the  modified  hollow  cylinders  appeared  to  give  the  best  overall 
increases  in  both  pitch  and  yaw  simulation.  An  increase  of  approximately  9  deg  in  upwash 
was  obtained  and  was  fairly  uniform  over  a  region  typical  of  an  inlet  capture  area.  By 
using  a  combination  of  flow  shaping  and  geometric  pitch,  the  capability  of  simulating 
20-deg  angle  of  attack  for  the  study  configuration  was  demonstrated.  The  yaw  simulation 
capability  was  also  demonstrated  with  the  modified  hollow  cylinders  by  using  a 
combination  of  flow  shaping  and  geometric  pitch. 

Although  the  pitch  simulation  obtained  was  5  deg  short  of  the  objective,  the  data 
obtained  show  that  this  testing  technique  can  significantly  increase  the  present  testing 
capability  of  the  AEDC  PWT-16T  for  testing  of  inlet/engine  systems,  particularly  those 
that  are  fuselage  side  mounted  and  are  located  forward  of  the  wing.  Further  study  is 
needed  before  it  can  be  assumed  that  the  technique  will  work  for  all  inlet/engine  locations. 
In  terms  of  modifications  or  additions  that  need  to  be  added  to  PWT-16T  to  utilize  this 
technique,  the  following  recommendations  are  made: 

1.  A  double  compressor  unit  equal  to  the  present  double  units  needs  to  be 
added  to  the  PWT-PES.  This  unit  is  required  to  cover  the  added  auxiliary 
flow  with  the  flow  shaping  devices  and  a  full-scale  inlet/engine  model. 

2.  A  test  cart  needs  to  be  modified  in  order  to  have  remote  control  of  the 
flow  shaping  devices  in  pitch,  yaw,  and  vertical  position  in  the  tunnel. 
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APPENDIXES 
I.  ILLUSTRATIONS 
II.  VORTEX-LATTICE  PROGRAM 
FOR  3  D  POTENTIAL  FLOW 
PROBLEMS 
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b.  Inlet  Performance  for  a  %-Scale  Model  (Ref.  1) 

Fig.  1  Typical  Performance  for  Highly  Maneuverable  Aircraft 
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Fig.  2  Aircraft  Study  Configuration  Installed  in  the  AEDC  16-ft  Transonic  Wind  Tunnel 
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Fig.  5  Mathematical  Model  of  Dual  Hollow  Cylinders  in  the  Wind  Tunnel 
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Fig.  7  Predicted  Sidewash  for  the  Dual  Hollow  Circular  Cylinders,  Mach  Number  0.9 
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Fig.  8  Mathematical  Model  of  the  Dual  Hollow  Circular  Cylinder  Configuration 
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Fig.  10  Mathematical  Model  of  the  Modified  Hollow  Cylinder  Configuration 
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Fig.  12  Mathematical  Model  of  the  Large  Turning  Vanes 
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Fig.  13  Local  Flow  Angles  for  Large  Turning  Vanes  in  Plane  of  Trailing  Edge 
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Fig.  16  General  Arrangement  of  Tunnel  IT  and  Supporting  Equipment 
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Fig.  17  Schematic  of  Tunnel  IT  Test  Leg 
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Orifice  Diameter  =  0.020  in. 
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Fig.  18  Flow  Angularity  Probe  and  Support  Mechanism 
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Fig.  19  Schematic 
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the  Dual  Hollow  Circular  Cylinder  Installation 
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Fig.  20  Dual  Hollow  Cylinders  installed  in  the  AEDC  1-ft  Transonic  Wind  Tunnel 
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•  Probe  Tip  Position  for  Strut  and  Roll  Mechanism  Zero 
Position  for  9.625  in.  Spacing  Outside  to  Outside 


Fig.  21  Probe  Survey  Grid  for  the  Dual  Hollow  Circular  Cylinder  Wind  Tunnel  Test 
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All  Dimensions  in  Inches 

•  Probe  Tip  Position  for  Strut  and  Roll  Mechanism  Zero 
Position  for  8.625  in.  Spacing  Outside  to  Outside 


Fig.  23  Probe  Survey  Grid  for  the  Modified  Hollow  Cylinder  Wind  Tunnel  Test 
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Fig.  24 


Schematic  of  the  Large  Turning  Vane  Installation 
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Fig.  25  Large  Turning  Vanes  Installed  in  the  AEDC  PWT-1T 
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•  Probe  Tip  Position  for  Strut  and  Roll  Mechanism  Zero 
Survey  Taken  at  Trailing  Edge  of  Vane  No.  2 

Fig.  26  Probe  Survey  Grid  for  the  Large  Turning  Vane  Wind  Tunnel  Test 
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Fig.  27  Schematic  of  the  Small  Turning  Vane  Installation 
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Fig.  28  Tunnel  Turning  Vanes  Installed  in  the  AEDC  PWT-1T 
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All  Dimensions  in  Inches 

•  Probe  Tip  Position  for  Strut  and  Roll  Mechanism  Zero 
Survey  Made  at  Trailing  Edge  of  Vanes 

Fig.  29  Probe  Survey  Grid  for  the  Small  Turning  Vane  Wind  Tunnel  Test 
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Mach  Number  =  0.9 
Experimental  Flow-Field  Data  - 
Theoretical  Flow-Field  Data  — 


Cylinder  Diameter  56  in.  Full  Scale, 
Cylinders  Spaced  154  in. 
Center-to-Center  Full  Scale 


Fig.  31  Experimental  and  Theoretical  Comparison  of  Local  Flow  Vectors  for  the  Dual  Hollow  Circular  Cylinders 
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Fig.  33  Experimental  and  Theoretical  Comparison  of  Local  Mach  Numbers  for  the  Dual 
Hollow  Circular  Cylinders 
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0,  deg 

Fig.  34  Experimental  Pressure  Coefficients  around  Cylinder  Number  1  of  the  Dual 
Hollow  Circular  Cylinder  Configuration 
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Cylinder  No.  1  -  a  ■  30  deg,  (3  -  0  deg 
Cylinder  No.  2  -  a  -  20  deg,  3  -  0  deg 


a.  Variation  in  Local  Mach  Number 


Survey  data  were 
not  obtained  for  2 
stations  at  Mach 
Number  0.6. 


b.  Change  in  Upwash 


Free-Stream  Mach  Number 


c.  Change  in  Sidewash 

Fig.  35  Variation  in  Experimental  Flow  Field  with  Free-Stream  Mach  Number 
for  the  Dual  Hollow  Circular  Cylinders 
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Mach  Number  -  0.9 

Experimental  Flow-Field  Data  — —  Cylinders  Spaced  138  in. 

Theoretical  Flow-Field  Data - -  Center-to-Center  Full  Scale 

0  =  0  510  20  30  deg 


Fig.  36  Experimental  and  Theoretical  Comparison  of  Local  Flow  Vectors  for  the  Modified  Hollow 
Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  20  deg) 
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Fig.  37  Experimental  and  Theoretical  Comparison  of  Upwash  and  Sidewash  Angles  for  the  Modified 
Hollow  Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  20  deg) 
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Cylinders  Spaced  138  in. 
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Fig.  38  Experimental  and  Theoretical  Comparison  of  Local  Mach  Numbers  for  the  Modified  Hollow 
Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  20  deg) 
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(Cylinder  No.  1) 


Fig.  39  Experimental  Pressure  Coefficients  around  Cylinder  No.  1  of  the  Modified  Hollow 
Cylinder  Configuration 
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a.  Variation  in  Local  Mach  Number 


Free-Stream  Mach  Number 

c.  Change  in  Sidewash 

Fig.  40  Variation  in  Experimental  Flow  Field  with  Free-Stream  Mach  Number  for  the  Modified 
Hollow  Cylinders  {Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  20  deg) 
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Mach  Number  =  0.9 

Experimental  Flow-Field  Data - —  Cylinders  Spaced  138  in. 

Theoretical  Flow-Field  Data - -  Center-to-Center  Full  Scale 

0  =  0  5 10  20  30  deg 


Fig.  41  Experimental  and  Theoretical  Comparison  of  Local  Flow  Vectors  for  the  Modified  Hollow 
Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  30  deg) 
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Mach  Number  -  0.9 
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Fig.  42  Experimental  and  Theoretical  Comparison  of  Upwash  and  Sidewash  Angles  for  the  Modified 
Hollow  Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  30  deg) 
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Fig.  43  Experimental  and  Theoretical  Comparison  of  Local  Mach  Numbers  for  the  Modified  Hollow 
Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  30  deg) 
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Free-Stream  Mach  Number 


a.  Variation  in  Local  Mach  Number 


b.  Change  in  Upwash 


c.  Change  in  Sidewash 

Fig.  44  Variation  in  Experimental  Flow  Field  with  Free-Stream  Mach  Number  for  the 

Modified  Hollow  Cylinders  (Cylinder  No.  1,  a  =  30  deg;  Cylinder  No.  2,  a  =  30  deg) 
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Fig.  46  Variation  in  Experimental  Local  Flow  Vectors  with  Yaw  Angles  for  Modified  Hollow  Cylinders 
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Mach  Number  -  0.9 
Experimental  Flow-field  Data  — — 
Theoretical  Flow-Field  Data - — 


8-  0  5 10  20  30  deg 
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Survey  Taken  at  Trailing  Edge  of  Vane  No.  2 

Turning  Vanes  Span  the  Tunnel,  96  in.  Apart  Full  Scale  (Fig.  24) 
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Fig.  46  Experimental  and  Theoretical  Comparison  of  Local  Flow  Vectors  for  Large 
Turning  Vanes 
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•  Experimental  Data 

Theoretical  Upwash,  e  =  Odeg 
(or  Sidewash  a  »  Odeg) 

(V) - Theoretical  Sidewash 

(or  Upwash) 


Survey  Taken  at  Trailing  Edge  of  Vane  No.  2 
Turning  Vanes  Span  the  Tunnel,  96  in.  Apart  Full  Scale 

Fig.  47  Experimental  and  Theoretical  Comparison  of  Upwash  and  Sidewash  Angies  for 
the  Large  Turning  Vanes 
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Fig.  48  Experimental  Surface  Pressure  Coefficients  for  the  Large  Turning  Vanes 
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Fig.  49  Experimental  and  Theoretical  Comparison  of  Local  Flow  Vectors  for 
Small  Turning  Vanes 
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Fig.  50  Experimental  and  Theoretical  Comparison  of  Upwash  and  Sidewash  Angles  for  the  Small  Turning  Vanes 


AEDC-TR-73-9 


Flow  Field  Survey 


-  High  Performance  Aircraft  Configuration 

(PWT-16T  Test  Data)  a  - 10  deg,  0  *  0  deg 

-  Dual  Hollow  Cylinders  (PWT-1T  Test  Data) 

Data  Adjusted  for  5-deg  Pitch  of  the  Inlet-Engine 
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Fig.  51  Comparison  of  Local  Flow  Vectors  in  the  AEDC  PWT-16T  with  the  Study  Configuration  at  an  Angle  of 
Attack  of  10  deg  with  Simulated  Local  Flow  Vectors  Taken  in  the  AEDC  PWT-1T 
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Flow  Field  Survey 

-  High-Performance  Aircraft  Configuration 

(PWT-16T  Test  Data)  a  -  19. 8  deg,  (5  ■  0  deg 

- Dual  Hollow  Cylinders  (PWT-1T  Test  Data) 

Data  Adjusted  for  12-deg  Pitch  of  the 
Inlet-Engine 


Position 


Fig.  52  Comparison  of  Local  Flow  Vectors  Taken  in  the  AEDC  PWT-16T  with  the  Study  Configuration  at  an 
Angle  of  Attack  of  19.8  deg  with  Simulated  Local  Flow  Vectors  Taken  in  the  AEDC  PWT-1T 
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Flow-Field 


- Experimental,  High-Performance  Aircraft  Configuration 

(PWT-16T  Test  Data) a  - 10 deg,  (3  -  4  deg,  Lee  Inlet 

- Theoretical,  Modified  Cylinders  Spaced  154  in.  Center-to-Center  Full  Scale 

(Data  Adjusted  for  10  deg  Geometric  Pitch  of  Inlet/Engine  System.) 
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Fig.  53  Comparison  of  Local  Flow  Vectors  Taken  in  the  AEDC  PWT-16T  with  the  Study  Configuration 
at  an  Angle  of  Attack  of  10  deg  and  an  Angle  of  Yaw  of  4  deg  with  Simulated  Local  Flow 
Vectors  Taken  in  the  AEDC  PWT-1T 
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_ (deq) —  Experimental,  High-Performance  Aircraft  Configuration 

<PWT-16T  Test  Data)  a  -  lOdeg,  P  -  4  deg,  Lee  Inlet 

(cleg) - Theoretical,  Modified  Cylinders,  154  in.  Center-to-Center  Full  Scale 

(Oata  Adjusted  for  lOdeg  Geometric  Pitch  of  Inlet/Engine  System) 
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Fig.  54  Experimental  and  Theoretical  Comparison  of  Upwash  and  Sidewash  Angles  for  the 
Study  Configuration  and  Modified  Hollow  Cylinders 
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Fig.  55  Comparison  of  Local  Flow  Vectors  Taken  in  AEDC  PWT-16T  with  the  Study  Configuration 
at  an  Angle  of  Attack  of  10  deg  and  an  Angle  of  Yaw  of  4  deg  with  Simulated  Local 
Flow  Vectors  Taken  in  the  AEDC  PWT-1T 
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Fig.  56  Installation  of  the  Inlet/Engine  Model  with  the  Dual  Hollow  Circular  Cylinders 
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Fig.  57  Installation  of  the  Inlet/Engine  Model  with  the  Modified  Hollow  Cylinders  in  the  AEDC  PWT-1T 
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o  Cylinder  No.  1,  a  ■  30  deg,  |3  ■  0  deg; 

Cylinder  No.  2,  a  -  20  deg,  |3  -  Odeg 
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Fig.  58  Increase  in  Upwash  and  Sidewash  Angles  due  to  the  Installation  of  the 
Inlet/Engine  with  the  Dual  Hollow  Circular  Cylinders 
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- Dual  Hollow  Cylinders 

- Modified  Hollow  Cylinders 


Fig.  59  Auxiliary  Weight  Flow  Required  to  Operate  the  Tunnel  with  Flow  Shaping 
Devices  and  Inlet/Engine  Installed 
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APPENDIX  II 

VORTEX-LATTICE  PROGRAM 
FOR  3-D  POTENTIAL  FLOW  PROBLEMS 

NOMENCLATURE 

a.  Speed  of  sound  at  critical  conditions 

C*  Defined  by  Eq.  (II-3) 

CP  Pressure  coefficient,  Eq.  (11-23) 

k  Ratio  of  specific  heats,  Cp/Cv 

M  Local  Mach  number,  Eq.  (11-22) 

M„  Free-strearn  Mach  number 

v*  Velocity 

t  Free-stream  velocity,  O  or  unit  velocity  in  X-direction 

U  { 

v  >  Velocity  components 

w  ) 

v^-rfk  Vector  dot  product  of  the  velocity  and  the  unit  normal  at  the  K'th  control 

point 

7  Strength  of  a  horseshoe  vortex 

17  Strength  of  a  horizontal  vortex  segment  of  the  model 

MNSUM  Total  number  of  independent  horseshoe  vortices  of  the  model 

GENERAL  INFORMATION 

The  program  is  written  in  FORTRAN  and  is  run  on  the  IBM  370/155. 

Purpose 

This  program  was  originally  written  at  Boeing  (Ref.  5).  At  ARO  it  was  modified 
several  times  by  several  different  people,  to  add  Goethert's  rule,  restart  capability,  etc. 
Use  of  the  program  became  complicated.  Some  of  the  options  were  no  longer  available 
but  still  were  included  in  the  input.  Internal  coding  was  patchwork  and  inefficient. 
Therefore,  the  program  has  been  completely  rewritten. 
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The  program  computes  velocities  and  streamlines  in  3-D  potential  flow  problems. 
The  boundary  conditions  to  be  specified  are  the  free-stream  velocity  and  the  vector  dot 
product  of  the  velocity  and  the  normal  to  the  boundary.  Compressibility  can  be  accounted 
for  using  Goethert's  rule. 

The  program  has  been  used  to  compute  subsonic  flow  around  wings  and ’airplanes 
and  to  compute  flow  around  various  configurations  in  wind  tunnels.  It  has  also  been  used 
to  compute  flow  fields  of  V/STOL  aircraft. 

Assumptions 

The  boundary  is  represented  by  a  lattice  network  of  vortices.  The  lattice  is  analyzed 
by  the  program  into  the  correct  number  of  independent  horseshoe  loops,  and  a  control 
point  is  located  on  the  boundary  for  each  loop.  Then,  the  strengths  of  the  loops  are 
calculated  so  that  the  boundary  conditions  are  satisfied  at  the  control  points.  When  the 
strengths  are  known,  the  velocity  at  any  point  in  the  field  can  be  computed  using  the 
Biot-Savart  law. 

Limitations 

The  boundary  is  represented  by  a  network  of  vortices.  The  velocity  is  very  large 
at  points  near  a  vortex.  Near  the  boundary  the  velocity  is  almost  wholly  the  effect  of 
one  or  two  vortices,  not  representative  of  the  real  boundary  in  that  locality;  therefore, 
the  velocity  is  not  accurate.  Moving  away  from  the  boundary,  the  velocity  is  the  effect 
of  several  vortices,  and  the  solution  becomes  accurate.  One  can  investigate  the  flow  closer 
to  the  boundary  by  making  a  finer  lattice;  however,  this  rapidly  increases  computer  time. 

The  program  is  for  potential  flow  problems;  however,  compressibility  can  be  accounted 
for  using  Goethert's  rule.  Good  results  have  been  obtained  for  Mach  numbers  up  to  0.9. 

The  program  is  limited  to  steady  flow. 

METHOD 

Introduction 

Consider  the  problem  of  calculating  the  subsonic  flow  about  a  wing.  To  use  this 
program,  the  wing  must  be  modeled  by  a  lattice  network  of  vortices  which  are  input 
to  the  program.  The  program  analyzes  the  lattice  into  the  correct  number  of  independent 
horseshoe  loops.  For  each  loop,  a  control  point  is  calculated  and  a  unit  vector  normal 
to  the  wing  at  the  control  point.  The  strengths  of  the  horseshoe  vortices  are  determined 
such  that  the  velocity  at  the  control  points  is  perpendicular  to  the  normal  vectors.  Once 
the  strengths  of  the  vortices  are  known,  the  velocity  at  any  point  can  be  computed  from 
the  Biot-Savart  law. 
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In  the  above  problem,  there  is  a  free-stream  velocity  and  the  boundary  condition 
that,  at  the  surface  of  the  wing,  the  flow  is  parallel  to  the  surface.  The  program  can 
also  be  used  for  problems  where  there  is  flow  through  a  surface  as  in  the  case  of  a 
fan-in-wing  configuration  (Ref.  6),  flow  in  a  wind  tunnel  with  porous  walls,  or  the  efflux 
tube  of  V/STOL  aircraft.  The  free-stream  velocity  can  be  set  to  zero  or  be  set  to  unit 
magnitude  in  the  direction  of  the  x-axis. 

The  Model 

A  model  is  input  to  the  program  in  several  parts,  divided  according  to  the  convenience 
or  experience  of  the  user.  A  part  is  illustrated  in  Fig.  II- 1.  The  figure  is  shown  flat, 
rectangular,  and  horizontal,  but  the  lattice  can  be  twisted,  stretched,  and  rotated  as  desired 
to  lie  on  the  boundary  surface.  The  number  of  points  and  the  coordinates  of  each  point 
are  input.  The  number  of  points  for  the  part  shown  would  be  input  as  follows: 

M  =  4  One  more  than  the  number  of  vertical  lines 

N  =  5  The  number  of  horizontal  lines 

ITR  =  3  The  number  of  points  along  a  horizontal  line 

minus  M,  called  extra  points 

For  any  part,  M  >  2,  N  3*  2,  and  ITR  >  0. 

The  coordinates  are  input  in  the  following  order.  First  for  J  =  1,  the  coordinates 
at  I  =  1,  2,  ...,  M  +  ITR  are  input.  Next  for  J  =  2,  the  coordinate  at  I  =  1  to  I  = 

M  +  ITR  are  input.  Then  the  same  for  J  =  3  and  so  on  to  J  =  N.  All  the  coordinates 

for  the  first  part  are  input,  then  the  second,  and  so  on  for  all  parts. 

The  part  shown  in  Fig.  II- 1  is  analyzed  by  the  program  into  12  horseshoe  vortices 
as  shown  in  Fig.  II-2.  Note  that  the  number  of  independent  horseshoe  vortices  is  (M 
-  1)  (N  -  1).  All  the  horseshoe  vortices  for  a  given  part  trail  off  from  the  last  point 
to  infinity  at  a  given  direction.  The  direction  is  input  for  each  part  via  the  variables  AX, 
AY,  and  AZ,  that  is,  the  angle  in  degrees  that  the  trails  make  with  the  X,  Y,  and  Z 
axes,  respectively.  A  trail  can  be  made  to  fall  along  a  curve  by  using  extra  points,  that 
is  ITR  >  0.  A  lattice  with  no  trailing  vortices  can  be  made  by  using  extra  points  to 

make  the  last  points  of  all  the  horizontal  lines  coincide;  thus,  all  the  trails  will  fall  along 

each  other  and  cancel  out.  That  is,  the  net  strength  of  the  trail  (singular  since  all  the 
trails  coincide)  will  be  zero. 

The  program  locates  a  control  point  in  the  center  of  the  first  4  points  of  each 
horseshoe  vortex,  represented  by  X's  in  Fig.  1 1-2  and  computes  a  unit  normal  at  those 
points.  If  there  is  flow  through  the  body,  then  the  vector  dot  product  of  the  velocity 
and  the  unit  normal  at  each  control  point  is  input.  They  are  input  in  the  following  order: 
First  for  J  =  1,  the  velocity  dot  normals  at  control  points  I  =  1,2,  ...,  M  -  1  are  input. 
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Next  for  J  =  2,  they  are  input  for  I  =  1  to  I  =  M  -  1.  Then  the  same  for  J  =  3  and 
so  on  to  J  =  N  =  1.  All  the  velocity  dot  normals  for  the  first  part  are  input,  then  the 
second,  and  so  on  for  all  parts. 

The  description  in  this  section  uses  terms  such  as  vertical,  horizontal,  rectangle,  etc., 
referring  to  the  figures.  It  is  to  be  remembered  that  an  actual  part  can  have  any  orientation, 
for  example,  rotated  90  deg,  rotated  180  deg,  stretched,  twisted,  and  have  a  completely 
different  appearance  than  the  figures,  determined  completely  by  the  coordinates  input. 

Theory 

The  velocity  at  a  point  induced  by  a  horseshoe  vortex  is  given  by  the  Biot-Savart 
law 


V  =  X/. 
4  V  r 


'  X  df 


(11-1) 


where  the  line  integral  is  over  the  vortex  filament,  r*  is  the  vector  from  the  point  to 
the  vortex,  and  ds  is  of  arc  length  magnitude  and  is  tangent  to  the  vortex  filament. 


The  velocity  at  any  point  is  equal  to  the  free-stream  velocity  plus  the  sum  of  the 
velocity  induced  at  that  point  by  each  horseshoe  vortex.  Therefore,  the  velocity  at  the 
K'th  point  is 


MM  SUM 
£ 
j=l 


(H-2) 


where  "vkj  is  the  velocity  induced  at  the  K'th  point  by  the  J'th  horseshoe  vortex  and 
MNSUM  is  the  number  of  vortices.  The  velocity  induced  by  a  vortex  is  directly  proportional 
to  its  strength,  so  define 


c4  • 

y 


(H-3) 


then  Eq.  (II-2)  becomes 


MNSUM 


v.  =  v„  + 


1  C..  y. 

;-i  ^3  y) 


3= 1 


(II-4) 


Consider  the  case  where  the  K'th  point  refers  to  the  K'th  control  point  and  take  the 
dot  product  of  both  sides  of  Eq,  (II-4)  with  Ttk  and  rearrange 


MNSUM 


.1  [q/  •  nk]yj  =  vk  .  nk 


-  v  •  n, 


(H-5) 


85 


AEDC-TR-73-9 


The  Vk'nt  is  either  0  or  input.  The  free-stream  velocity  is  either  0  or  input.  The  n^ 
are  computed.  (*kj  is  the  velocity  induced  at  the  K’th  control  point  by  the  J’th  vortex 
when  its  strength  is  one;  therefore,  it  can  be  computed.  Thus,  everything  in  Eq.  (11-5), 
except  7j,  is  either  known  or  can  be  calculated,  and  Eq.  (11-5)  is  a  linear  algebraic  set 
of  equations  which  can  be  solved  for  the  strengths  of  the  vortices. 


When  strengths  are  known,  the  velocity  at  any  point  K  can  be 
(II-2).  Streamlines  are  calculated  from  velocities  by  integrating 

obtained  from  Eq. 

-*  -4 

dr  v 

dr  v 

(11-6) 

Symmetry 

Consider  a  system  that  is  symmetrical  with  respect  to  a  plane.  Number  half  of  its 
vortices  from  1  to  MNSUM  and  the  other  half  from  MNSUM+1  to  2*MNSUM  where 
2*MNSUM  corresponds  to  1,  2*MNSUM-I  corresponds  to  2,  and  so  on  to  MNSUM+1 
corresponds  to  MNSUM.  Equation  (1 1-5)  becomes 

2*  MNSUM  _  _ 

[Ckj  •  "k^j  =  vk  *  "k  “  v~  •  "k 

j=l 

(11-/) 

Define 

i  =  2*MNSEM+l-j 

(11-8) 

l  =  2*\1NSUM+1-K 

(11-9) 

By  symmetry, 

it  is  seen  that  v^  is  parallel  to  the  plane  of  symmetry  and  that 

I 

II 

X 

(11-10) 

vf  •  =  vk  •  ”k 

(11-11) 

v  *  n  p  =  v  •  n , 

oo  E  oo  ^ 

(11-12) 

— *  — > 

CU  •  "f  “  Ckj  •  "k 

(11-13) 

By  changing  the  free  index  to  m,  combining  corresponding  terms  on  the  left  hand  side, 
and  using  Eqs.  (11-8)  and  (11-10),  then  Eq.  (11-7)  becomes 


MNSl’M  _ 

.Jj  [<Cnj  -  Cmi)  •  njyj  =  vm  •  .  nm  (11-14) 
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It  is  seen  that  Eq.  (II- 1 4)  yields  the  same  equation  for  m  =  k  and  for  m  =  C  Therefore, 
one  can  solve  Eq.  (11-14)  for  7,  for  j  from  I  to  MNSUM  and  use  Eq.  (II- 10)  for  the 
rest  of  the  7's. 

Computing  the  coefficients  for  a  symmetrical  system  takes  only  half  as  long  since 
the  coefficients  for  only  half  of  the  equations  need  he  computed.  Solving  the  system 
takes  only  about  one-eighth  as  long,  since  the  time  to  solve  the  equations  goes  up 
approximately  like  the  cube  of  the  number  of  equations.  The  program  can  be  used  with 
no  symmetry,  symmetry  with  respect  to  the  XY  plane,  symmetry  with  respect  to  the 
XZ  plane,  or  symmetry  with  respect  to  both  the  XY  and  XZ  plane.  The  latter  ease  requires 
only  one-fourth  as  long  to  compute  the  coefficients  and  about  1  /64th  as  long  to  solve 
the  equations  as  an  equivalent  nonsymmetrical  case. 

Goethert's  Rule 


All  input  and  output  to  the  program  is  in  the  physical  reference  system.  Internally, 
however,  a  transformation  is  made  to  account  for  compressibility  according  to 


v'l  -  Mr. 


(11-15) 


(11-16) 

(11-17) 


Velocity  components  in  the  transformed  plane  (u\  v',  w') 
transformed  back  to  the  physical  plane  according  to 


u '  -  M  - 


1  -  M2 


y/T-tf 

V  OO 


\/l  -  M  2 

*  OO 


are  obtained  and  then  are 

(11-18) 

<U-1  9) 

(11-20) 


Transformations  in  the  x-direction  are  different  than  transformations  in  the  y  and  z 
directions,  because  v„  is  in  the  x  direction.  These  transformations  are  according  to 
Goethert's  rule  as  given  in  Ref.  7. 


Equations  for  Local  Flow  Data  That  Are  Output  by  the  Program 

The  flow  data  output  by  the  program  for  a  given  point  are  x.  y,  z,  u,  v,  w ,  |v|, 
M,  tan-1  v/u,  tan'1  w/u,  CP,  and  M  -  M„.  The  coordinates  of  the  point  are  x,  y,  z,  and 
the  components  of  the  velocity  are  u,  v,  w.  The  magnitude  of  the  velocity  is 
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|v|  =  \/u“  +  V2  +  w2 

The  local  Mach  number  is  computed  by 


(11-21) 


(11-22) 


The  pressure  coefficient  by 

CP  =  1  -  |  v|  2  (11-23) 

The  How  angularity  is  determined  by  tan-1  v/u  and  tan-1  w/u.  These  parameters  are 
computed  as  indicated  as  is  M  -  M^. 

Equation  (11-22)  is  derived  from  Eqs.  (50)  and  (58)  of  Ref.  8.  Since  the  free-stream 
velocity  is  of  unit  magnitude,  the  above  mentioned  Eq.  (50)  becomes  in  the  free-stream 


Equation  (58)  of  Ref.  8  is 

_J_  hi2 

9  k  +  1  a2 

M2  =  - ^ -  (11-25) 

k  -  1  |  v|  2 

1  + 

k  +  1  2 

a* 

By  eliminating  at  and  solving  for  M,  one  obtains  Eq.  (11-22). 

Effect  of  a  Vortex  Segment  or  a  Vortex  Ray 

Once  7's  have  been  obtained  it  would  be  very  inefficient  to  compute  Eq.  ( 1 1- 1 )  for 
each  horseshoe  vortex  and  add  the  results  to  calculate  a  velocity.  Instead,  the  net  strength 
of  each  segment  and  ray  is  computed,  and  Eq.  (11-1)  is  evaluated  just  once  over  each 
segment  and  ray  of  the  model.  Note  that  the  strength  of  each  vertical  segment  (Fig.  11-2) 
is  the  strength  of  its  respective  horseshoe  vortex.  The  strength  of  a  horizontal  segment, 
however,  is  the  algebraic  sum  of  the  strengths  of  each  horseshoe  vortex,  of  which  the 
segment  is  a  part.  Let  n; j  be  the  net  strength  of  a  horizontal  horseshoe  segment  where 
j  denotes  the  horizontal  line  and  i  which  segment  of  the  line.  It  is  seen  from  Fig.  1 1-2 
that 


z)  j  j  =  £l  +  £2  +  (=3 


(Il-26a) 


88 


AEDC-TR-73-9 


where 


Cl  = 
C2  = 
C3  = 


•  o 

I  0 

r« 


1° 


i  >  1 
i  =  1 

j  >  1 
j  =  1 
j  <  N 
j  =  N 


(Il-26b) 


(ll-26c) 


(11-26(1) 


Since  the  model  is  made  up  completely  of  vortex  segments  and  vortex  rays  (the 
trailing  vortices),  only  two  formulas  are  needed  by  the  program  to  evaluate  Eq.  (II-l), 
one  where  the  curve  C  corresponds  to  a  line  segment  and  the  other  when  it  is  a  ray. 
In  Fig.  11-3,  the  line  segment  extends  from  point  A  to  point  B,  and  C  is  the  point  at 
which  the  velocity  is  desired.  Equation  II- 1  results  in 

v  =  (cost),  +  cos  (?9)u  (11-27) 

4?7a  1  L 

where  "u  is  a  unit  vector,  which  for  the  case  shown  in  Fig.  11-3,  is  directed  into  the  paper. 
The  formula  for  a  ray  is  easily  obtained  by  letting  Oj  approach  zero.  Thus,  it  is  identical 
to  Eq.  (11-27)  except  cos  is  replaced  by  1. 

RESULTS 

Figure  14  shows  the  model  that  was  used  in  the  sample  problem.  The  model  is 
symmetrical  with  respect  to  the  XY  plane,  and  only  one-half  of  the  model  is  shown. 
Figure  15  contains  contour  plots  of  the  flow  angularity  computed  by  the  program  for 
this  model.  Printed  output  for  this  shot  follows  the  program  listing. 

INPUT  PREPARATION  AND  OUTPUT  DESCRIPTION 
Card  Input  Data 

Card  input  to  the  program  is  indicated  below  in  "pseudo-FORTRAN."  The  format 
of  each  card  is  given  in  parenthesis  after  READ.  Definitions  of  variables  follow  the 
"FORTRAN." 

1  INPUT  TO  VORTEX  LATTICE  PROGRAM 

2 

3  READ  (213)  I  REST,  MINGO 

4  READ  (A80)  TITLE 

5  READ  (413)  NW,  LXY,  LXZ,  L900 

6  READ  (El 2.0)  RER 
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7  RHAD  (2hl2.Cn  M„,  k 

8 

9  DO  1 0  K  =  I ,  NVV 

10  READ  (313.  3X,  3HI2.0),  M,  N,  ITR,  AX.  AY,  AZ 
I  1 

12  BODY  COORDINATES 

13  DO  16  K=  I  ,NW 

14  DO  16  J=l  ,N(K) 

15  DO  16  1  =  I.M(K)+ITR(K) 

16  RHAD  (3H12.0)  XB,  YB,  ZB 

17 

18  VELOC  ITY  DOT  NORMALS 

10  IF  (LO00-0)  GO  70  23 

20  READ  (6EI2.0)  «("vrf.  I=I,M(K)-I),  J=I.N(K)),  K=1.NW) 

21 

22  VELOCITIES 

23  READ  (13)  NV 

24  IF  (NV=0)  GO  TO  29 

25  DO  26  J=  |,NV 

26  READ  (3EI2.0)  XV,YV,ZV 

27 

28  STREAMLINES 

29  READ  (LONS 

30  IF  (NS=0)  GO  TO  35 

31  DO  33  J=  l.NS 

32  READ  (6E  12.0,18)  XO,YO,ZO,DSO,AN,AX,NX 

33  READ  (6E12.0)  XN,XX,YN.YX,ZN,ZX 

34 

35  END  OF  CARD  INPUT 


Definitions  of  above  variables  follow: 

I  REST  0  Initial  shot 

1  Restart  in  SETM 

2  Restart  in  LSYSEQ 

3  7's  are  known 

MINGO  Time  in  minutes  the  program  is  to  run. 

NW  Number  of  model  parts 

LXY  \  Symmetry  indicators 

LXZ  )  0  indicates  no  symmetry 

I  indicates  symmetry 
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L900  0  _vao=u,  where  u  is  unit  velocity  in  x-direction,  ”v'Ti=0 

1  v<XJ=  u ,  vn  are  read 

2  "v^O,  rii  are  read 

RER  If  a  point  is  within  tlris  distance  of  a  line  determined  by  a 

vortex  segment,  the  effect  of  that  segment  is  set  to  zero. 
RF.R>0. 

M„  Free-stream  Mach  number 

k  Ratio  of  specific  heats.  If  left  zero,  it  defaults  to  1.4. 

M  One  more  than  the  number  of  vertical  lines.  \\>2. 

N  The  number  of  horizontal  lines.  N>2. 

ITR  The  number  of  extra  points  -  i.e,,  the  number  of  points  along 

a  horizontal  line  minus  M.  ITR>0. 

AX')  Angle  in  degrees  the  trailing  vortices 

AY  ^  of  the  part  make  with  the 

AZ  *  X,  Y,  and  Z  axis,  respectively. 

XB  ) 

YB  ^  Coordinates  of  points  on  the  model 

ZB  > 

vm  If  L.900^0,  then  the  "vii  must  be  read  in,  6  to  a  card. 

NV  Number  of  points  at  which  velocities  are  to  be  calculated 

XV  I  Coordinates  of  the  points  at  which  velocities  are  to  be 

YV  >  calculated 

ZV  ' 

NS  Number  of  streamlines  to  be  calculated 

XO  | 

YO  Coordinates  of  starting  point  of  streamline 

ZO  ‘ 

DSO  Absolute  value  of  DSO  is  maximum  step  size  used  in  calculating 

streamline.  It  is  also  the  length  of  the  initial  step  size. 

The  sign,  plus  or  minus,  determines  whether  the  downstream 
or  upstream  streamline  is  to  be  calculated. 
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AN 

AX 

NX 

XN,XX  | 
YN,YX  r 
ZN,ZX  ) 

Read  and 

Unit  5 
Unit  6 

Unit  10 

Unit  1  1 

Unit  12 


Minimum  angle  in  degrees.  If  the  angle  between  velocities 
at  the  beginning  of  a  step  and  at  the  end  of  a  step  is  less 
than  AN,  then  the  step  size  is  increased  unless  it  is  already 
at  its  maximum  value  which  is  the  absolute  value  of  DSO. 

Maximum  angle  in  degrees.  If  the  angle  between  velocities  at 
the  beginning  of  a  step  and  at  the  end  of  a  step  exceed  AX,  then 
the  step  size  is  decreased. 

Maximum  number  of  steps.  Calculation  of  the  streamline  terminates 
after  this  number  of  steps. 

Minimum  value  and  maximum  value  for  X,Y,  and  Z.  If  the 
streamline  extends  past  one  of  these  values,  calculation  of 
the  streamline  terminates. 

Write  Units 

Card  Reader  -  See  previous  section  for  card  input  data. 

Printer  -  All  output  is  labeled  using  the  same  nomenclature 
as  this  report. 

Direct  Access  Unit  -  This  unit  is  used  only  in  LSYSEQ.  If 
it  is  known  that  a  shot  will  terminate  in  SETM  because  of 
time,  then  Unit  10  need  not  be  defined. 

System  of  Equations  -  It  is  written  in  SETM  and  read  by  LSYSEQ. 
It  IREST>2,  then  Unit  1 1  need  not  be  defined.  Each  logical 
record  of  Unit  1 1  consists  of  one  row  of  coefficients  of  the 
system  of  equations  and  the  right  hand  constant.  SETM  computes 
a  row,  then  writes  it  on  11;  computes  another  row  and  writes 
it,  and  so  on.  After  the  last  record,  that  is  either  all 
coefficients  have  been  computed  or  else  SETM  must  terminate 
because  of  time,  in  either  case,  Unit  I  1  is  endfiled.  See 
Unit  21  for  more  information. 

When  LSYSEQ  terminates  because  of  time,  it  writes  the  partially 
solved  system  of  equations  out  on  Unit  12.  When  the  solution 
is  complete  LSYSEQ  rewinds  12,  writes  without  format  j=l  to 
MNSUM,  and  then  endfiles  12.  Note  that  either  one  or  the 
other,  either  the  partially  solved  system  of  equations  or  else 
the  solution  is  written  on  12.  Not  both.  The  7's  are  read 
from  Unit  12  by  AIRFL2. 
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Unit  13  Velocities  -  If  it  is  known  that  a  shot  will  terminate  in  SETM 

or  LSYSF.Q  because  of  time,  then  Unit  13  need  not  be  defined. 
Otherwise,  Unit  13  must  be  defined,  whether  or  not  any 
velocities  are  requested.  USERV  computes  the  velocity  at  a 
point  and  writes  a  record  on  Unit  13;  computes  the  velocity 
at  the  next  point  and  writes  a  record  on  13.  and  so  on  until 
all  velocities  requested  have  been  computed.  It  then  writes 
out  a  record  of  -l's  and  then  endfiles  13  even  if  no  other 
records  were  written,  that  is  even  if  no  velocities  were 
requested.  Each  record  is  written  without  format  thus: 

XV.  YV,  UV,  U,  V,  W,  (v|,  M,  tail"1  (v/u).  tan-'tw/u). 

CP,  Nl-M^.  These  are  the  same  velocity  data  that  are  printed. 

The  program  never  reads  Unit  13.  The  velocity  data  are 
saved  on  Unit  13  for  the  convenience  of  the  user  in  case 
further  processing  is  required,  for  example  plotting. 

Unit  14  Streamlines  -  If  it  is  known  that  a  shot  will  terminate  in 

SETM  or  LSYSEQ  because  of  time,  then  Unit  14  need  not  be 
defined.  Otherwise,  it  must  be  defined,  even  if  no  streamlines 
are  requested.  STREAM,  called  by  USERS,  writes  one 
record  on  Unit  14  for  each  point  along  a  streamline.  This 
record  is  identical  to  that  written  lor  velocities,  that 
is  each  record  is  written  without  a  format  thus:  XS.  YS,  ZS. 

U,  V.  W,  (vl,  M,  tan1  (v/u),  tan1  (w'u),  CP.  M-M„  where 
XS,  YS,  ZS  are  the  coordinates  of  a  point  on  the  streamline. 

After  computation  of  a  streamline  is  finished,  a  record 
of  -l's  is  written.  After  all  streamlines  requested  have 
been  completed,  Unit  14  is  endfiled.  It  is  endfiled  even  if 
no  records  were  written,  that  is  even  if  no  streamlines  were 
requested.  The  program  never  reads  Unit  14.  Streamline  data 
are  saved  on  Unit  14  for  the  convenience  of  the  user,  in 
case  further  processing  is  required,  for  example  plotting. 

Unit  21  When  restarting  in  SETM  or  in  USYSEQ,  the  tape  written  on  Unit 

11  in  the  previous  shot  must  be  mounted  on  Unit  21.  Otherwise. 
Unit  21  need  not  be  defined.  When  restarting  in  SETM,  SETM 
reads  a  record  from  Unit  21  and  writes  it  on  Unit  11:  reads 
the  next  record  and  writes  it,  and  so  on  until  it  reaches  end 
of  file  on  Unit  21.  It  then  proceeds  computing  the  coefficients 
a  row  at  a  time  and  writing  them  on  Unit  1  1.  When  restarting 
in  USYSEQ,  before  LSYSEQ  is  called.  SETM  reads  the  system  of 
equations  from  Unit  21  and  writes  then  on  Unit  11. 

Unit  22  When  restarting  in  LSYSEQ,  the  tape  written  on  Unit  12  in  the 

previous  shot  must  be  mounted  on  Unit  22.  Otherwise.  Unit  22 
need  not  be  defined. 
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OPERATING  INFORMATION 
Restart  Capability 

The  time  in  minutes  a  shot  is  to  run  is  input  via  MINGO.  At  the  very  start,  the 
program  calls  an  assembler  language  subroutine,  TIMEX.  TIMEX  sets  the  FORTRAN 
variable,  JRUB  equal  to  zero,  and  at  the  end  of  the  allotted  time  sets  JRUB  equal  to 
one.  The  variable  JRUB  is  checked  by  the  subroutines  SETM  and  LSYSEQ  at  the  end 
of  each  row  of  computation.  If  time  has  run  out,  they  terminate  and  can  be  restarted. 
Note  that  the  routines  check  at  the  end  of  a  row  of  computation,  then  must  write  restart 
information  on  tape,  thus  the  shot  will  run  a  little  longer  than  the  time  input  MINGO. 
This  must  be  considered  because,  if  the  program  is  killed  by  the  operating  system  instead 
of  terminating  normally,  then  it  might  be  impossible  to  restart  it.  Also,  note  that  the 
program  checks  only  the  time  while  calculating  the  coefficient  matrix  (SETM)  and  while 
solving  the  system  of  equations  (LSYSEQ).  Once  the  system  of  equations  have  been  solved 
(7's  are  obtained),  the  program  will  start  calculating  the  velocities  and  streamlines  ordered. 
It  will  continue  until  killed  or  until  finished.  However,  once  7's  are  obtained,  one  can 
always  restart,  calculating  velocities  at  left  over  points  or  restarting  a  streamline  where 
it  left  off. 

Running  Time 

The  time  in  minutes  for  a  run  can  be  roughly  estimated  as  follows: 

Time  in  SETM  %  Cl*  ( KS*MNSUM)**2 
Time  in  LSYSEQ  *  C2*  MNSUM**3 
Time  in  USERV  =  C3*  KS*  MNSUM  *NV 
Time  in  USERS  *  2*  C3*KS*MNSUM*NS*NX 


where 


1  no  symmetry 

KS=  2  symmetry  with  respect  to  one  plane 
4  symmetry  with  respect  to  two  planes 

NX  is  the  maximum  number  of  steps  for  a  streamline  and  NX  is  the  average  of 
the  NX's  for  all  streamlines. 

For  the  370/155,  the  constants  are 

Cl  =  1.8  X  10-4 
C2  =  1.7  X  10-6 
C3  =  2.6  X  10-5 
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Also  useful  in  estimating  time  is  the  fact  that  SETM  requires  approximately  the  same 
amount  of  time  for  each  row  of  computation.  However,  LSYSEQ  requires  less  and  less 
time  for  successive  rows. 

End  of  Shot  Indication 

The  last  line  of  printed  output  gives  the  reason  for  termination. 

For  example: 


1.  371  OF  430  IN  SETM 

2.  106  OF  430  IN  LSYSEQ 

3.  JOB  COMPLETED 

4.  SINGULAR  MATRIX 

5.  DIMENSIONED  BY  11000,  REQUIRES  12000 

6.  ZERO  NORMAL  3  7  2 

In  Example  1,  the  program  should  be  restarted  with  IREST=1.  In  Example  2,  the 
program  should  be  restarted  with  IREST=2.  In  Example  3,  all  requested  work  was 
completed.  If  additional  velocities  or  streamlines  are  desired,  one  can  restart  with  1REST=3. 
The  only  times  Example  4  has  occurred  were  for  bad  or  wrong  modeling.  Example  5 
can  be  eliminated  either  by  reducing  the  model  or  increasing  the  dimension  (see  the  next 
section  for  redimensioning  information.)  Best  use  of  storage  is  made  when  all  model  parts 
are  approximately  the  same  size,  that  is  all  parts  have  approximately  the  same  M.  N, 
and  1TR.  Thus,  Example  5  can  sometimes  be  eliminated  by  a  more  judicious  division 
of  the  parts.  Example  6  is  caused  by  a  modeling  error  where  the  first  four  points  of 
a  horseshoe  vortex  are  colinear.  The  three  integers  following  the  diagnostic  locate  which 
horseshoe  vortex  is  in  error. 

PROGRAMMING  INFORMATION 

Flow  Diagram  and  Modification  Hints 

A  flow  diagram  of  the  program  which  shows  where  each  subroutine  is  called  is  shown 
in  Fig.  11-4.  The  diagram  goes  from  left  to  right  and  from  top  to  bottom.  For  example, 
USERV  calls  VELLAB,  then  VELOCY,  then  VELOUT.  but  does  not  call  STREAM.  USERS 
calls  VELLAB,  then  STREAM,  and  STREAM  calls  VELOCY,  then  VELOUT. 

The  program  was  written  modular  in  form  for  easy  modification.  Figure  11-4  and 
the  following  section  should  be  helpful  if  modification  is  necessary.  As  an  example  of 
a  modification  which  has  been  made:  The  coordinates  for  one  large  model  were  calculated 
by  another  program  and  written  on  tape.  BODY  1  was  then  modified  to  read  the  coordinates 
off  the  tape  instead  of  cards. 
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Another  example:  It  was  easier  to  program  which  points  velocities  were  desired  rather 
than  prepare  a  deek  with  the  coordinates.  Therefore,  USERV  was  modified  by  replacing 
the  READ  statement  with  statements  producing  the  coordinates. 

Another  easy  change  which  could  be  made:  Suppose  some  other  parameters  which 
were  functions  of  the  components  u,  v,  and  w  were  desired  instead  of  or  in  addition 
to  the  parameters  now  calculated,  M,  |vj,  M-M„,  etc.  Then  one  need  only  modify  VELOUT 
to  compute  and  output  the  desired  quantities  and  modify  VEL.LAB  to  label  them  properly. 

To  redimension  the  program,  only  two  statements  need  be  changed.  They  are  the 
dimension  statement  in  the  main  program  and  the  integer  statement  in  subroutine  SCRIMP. 
SCRIMP  computes  how  much  storage  is  required  by  the  model  then  checks  to  see  if 
the  storage  is  available.  This  is  the  purpose  of  the  integer  statement  in  SCRIMP. 

Brief  Description  of  the  Subroutine 

TIMEX  is  an  assembler  language  routine.  When  called,  it  sets  the  FORTRAN  variable 
JRUB  equal  to  zero,  sets  a  clock,  then  returns.  Execution  of  the  program  continues  as 
usual.  When  the  allotted  time.  NSEC,  has  elapsed,  then  the  clock  sets  JRUB  equal  to 
one.  The  variable  JRUB  is  checked  by  SETM  and  LSYSEQ  after  each  row  of  computation. 
If  zero,  they  continue.  If  one,  they  terminate. 

TIMEX  must  be  supplied  by  the  user  according  to  the  particular  timing  technique 
of  the  computer  used.  Its  function  is  simple  and  clear.  It  sets  JRUB  to  zero  when  called, 
and  JRUB  must  change  to  one  after  the  elapsed  time.  NSEC.  One  could  if  necessary, 
by  sacrificing  the  restart  capability,  just  write  a  FORTRAN  subroutine  which  sets  JRUB 
to  zero  and  just  leave  it  zero.  Alternately,  one  could  modify  SETM  and  LSYSEQ  so  that 
instead  of  checking  JRUB,  they  compute  a  given  number  of  rows,  then  terminate. 

SCRIMP  computes  how  much  storage  is  needed  by  each  array.  The  storage  for  each 
array  is  allocated  in  the  single  array  X  of  the  main  program.  If  the  storage  requirement 
exceeds  the  dimension,  then  SCRIMP  writes  a  message  and  terminates  the  shot. 

CM  AIN  is  called  by  the  main  program  with  all  arguments  being  the  X  array  subscripted 
with  the  beginning  location  of  the  corresponding  dummy  argument.  These  beginning 
locations  were  supplied  by  SCRIMP.  CHAIN,  in  turn,  calls  the  subroutines  as  shown  in 
Fig.  1 1-4  with  the  proper  arguments. 

BODY  1  reads  the  body  coordinates,  sets^  according  to  L900,  and  either  reads  vm 
or  sets  them  to  zero  according  to  L900. 

A 1 R F L I  prints  the  input  and  sets  the  value  of  some  constants.  It  also  replaces  AX, 
AY,  and  AZ  with  their  cosines  producing  a  unit  vector  in  the  direction  of  the  trailing 
vortices. 
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STRECH  stretches  the  model  according  to  Goethert's  rule. 

BODY 2  computes  and  outputs  the  control  points  and  unit  normal  vectors. 

SETM  computes  the  coefficients  and  right  hand  constants  of  the  system  of  equations 
represented  by  Eq.  (11-5).  This  is  done  one  row  at  a  time  and  is  output  on  Unit  11. 

AIRFL2  calls  LSYSEQ  to  effect  the  solution  of  the  system  of  equations.  If  the  system 
is  singular  then  AIRFL2  writes  a  message  and  terminates.  Otherwise.  A1RFL2  reads  the 
7's  from  Unit  12  and  computes  the  tj’s.  It  prints  the  7's  and  p's  and  returns.  If  1REST>2. 
then  it  skips  calling  LSYSEQ  and  begins  by  reading  the  7's,  then  proceeds  as  before. 

USERV  reads  NV,  then  for  each  point:  (1)  reads  the  coordinate  of  the  point,  (2) 
obtains  the  velocity  at  the  point  by  calling  VELOCY.  then  (3)  outputs  by  calling  VELOUT. 
Before  returning  it  writes  a  record  of  -l's  on  Unit  13  then  endfiles  13. 

USERS  reads  NS.  then  for  each  streamline  requested:  (1)  reads  the  streamline  data, 
(2)  calls  stream,  then  (3)  writes  a  record  of  -l's  on  Unit  14.  Before  returning,  it  end- 
files  14. 

LSYSEQ  uses  direct  access  and  solves  a  large  system  of  equations.  If  direct  access 
is  not  available,  then  LSYSEQ  must  be  replaced  by  some  other  subroutine  using  a  different 
method.  Such  a  routine  is  usually  available  in  most  computer  centers.  The  system  of 
equations  is  on  Unit  1  I,  without  format,  with  one  row  of  coefficients  and  the  right  hand 
constant  per  record.  LSYSEQ  reads  Unit  11  and  solves  the  system  of  equations.  It  then 
writes  the  solution  on  Unit  12  in  one  record  without  format. 

DAF1 LE  is  an  assembler  language  routine  which  is  equivalent  to  the  DEFINE  FILE 
statement.  It  was  necessitated  by  the  requirement  that  the  arguments  of  the  DEFINE 
FILE  statement  be  constants.  DAF1LE  performs  exactly  the  same  function  as  DEFINE 
FILE  except  its  arguments  need  not  be  constants. 

STREAM  performs  a  numerical  integration  of  Eq.  (11-6)  and  outputs  via  VELOUT. 
The  method  used  is  improved  Euler. 

VELOCY  computes  the  velocity  at  a  point. 

VELOUT  computes  fvj,  M,  CP,  etc.,  and  prints.  It  also  outputs  on  a  designated  unit. 

VELLAB  is  an  entry  of  VELOUT  and  just  prints  a  label  for  the  printed  output 
of  VELOUT. 

SEGM  and  RAY  compute  the  velocity  induced  by  a  vortex  segment  and  a  vortex 
ray,  respectively. 
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VFYM  and  ZhMZ  replace  the  coordinates  of  the  model  by  the  coordinates  of  its 
reflection  with  respect  to  the  XZ  plane  and  the  XY  plane,  respectively. 


Fig.  11-1  Illustrative  Model  Cart 


98 


AE  DC-T  R-73-9 


Fig.  11-4  Flow  Diagram  of  Program 


Program  Listing  and  Sample  Run 

Following  is  a  program  listing  and  a  sample  run.  Figure  14  shows  the  model  of  the 
sample  run.  Figure  15  contains  contour  plots  of  the  flow  angularity. 
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G  LEVEL  20  MAIN  DATE  «  T3018  12/02/34 

C  SEP00211  DC  TODD  08/10/71  PW5146  RL  PALKD  B  1 

C  TIDY  VORTEX  LATTICE  PGM  B  2 

D 1MENSI0N  XI 11000)  B  3 

D (MENS  I  ON  Mil)  B  4 

EQUIVALENCE  (XII), Mill)  B  5 

COMMON  /CMMNA/  MI NGO , JRUB . I  RE S T , MNSUM .NWCODE , LXY , LX l , L900  B  6 

COMMON  /CMMNB/  RER , T I TLE I  20) , F S . AM, AG AM, AMS . F 1 . F2  B  7 

COMMON  /CMMNC  /  I  V  ,  J  V,  KV  ,L  TOT  ,  L  A  ,  LB  ,LC  ,L  D  ,LE  ,  LF  ,  LG  ,  LH  ,L  I  ,  L  J  ,L  K  ,L  L  ,  L  B  8 

1M,LN,L0,LP,LQ,LR,LS  B  9 

DIMENSION  LRAI1)  B  10 

EQUIVALENCE  I LR All), LA)  B  11 

READ  (5,40)  IREST, MINGO  B  12 

N  SEC®  69*M I NGO  B  13 

CALL  TIMEX  INSEC, JRUB)  B  14 

READ  (5,30)  TITLE  B  15 

READ  (5,40)  NWCODE ,LXY,LX2,L900  B  16 

READ  (5,50)  RER  ,F  S 

READ  (5,50)  AM, AGAM  B  18 

IF  (AGAM.EQ.O.)  AGAM«1.4  B  19 

LA=NWCDDE  B  20 

DO  10  K* 2,6  B  21 

10  LRAIK )«LRAtK-l)*NWCDDE  B  22 

DO  20  K«1 , NWCUDE  B  23 

20  READ  (5,60)  Ml K I , M( LA 4K ) , M( LB +K ) ,X ( LC«-K ) , X ( LD+K ) , X ( L E+K )  B  24 

CALL  SCRIMP  (Ml  1) ,M(LA+1)  ,M(LB+1) )  B  25 

CALL  CHAIN  I M I  1 ) , MI  LA ) ,M( LB  I , X( LC ) , X  I LD ) , X  I L E ) , XI LF ) , X I  LG ) , X I LH ) , X  B  26 

11  LI ),XILJ) ,X(LK) ,X(LL) ,X( LM)  ,X(LN) ,X(LO)  ,X (LP),X (LU) ,X(LR) )  B  2  7 

WRITE  (6,70)  B  28 

STOP  B  29 

30  FDRMAT  (20A4)  B  31 

40  FORMAT  (2413)  B  32 

50  FORMAT  (6E12.0)  B  33 

60  FDRMAT  I  31 3 , 3X , 3E 12 . 0)  B  34 

70  FORMAT  I 14HAJDB  COMPLETED)  B  35 

END  B  36- 


G  LEVEL  20  SCRIMP  DATE  »  7301B  12/02/34 

SUBROUTINE  SCRIMP  (M.N.ITR)  C  1 

DIMENSION  Mil),  Nil),  I TR ( 1 )  C  2 

COMMDN  /CMMNA/  KL, JRUB, (REST, MNSUM, NWCODE  C  3 

COMMDN  /CMMNC/  I  V, J V, K V ,L TOT , LI  1 )  .  C  4 

INTEGER  NDIM/UOOO/ 

MNSUM«0  C  5 

IV»0  C  6 

JV»0  C  7 

K  V*  0  C  8 

DO  10  K«1,NWCDDE  C  9 

MNSUM*MNSUM4(MIK)-1)*IN(K)-1I  C  10 

(VMAXOI  I  V.MIlOUTRlKI  )  C  11 

JV=MAXO( JV.MIK) )  C  12 

K  V»H  A  XOI  K  V,  N(  K I  )  C  13 

10  CONTINUE  C  14 

DO  20  K»l,6  C  |5 

20  L(K)®LIK)+1  C  16 

LA»( V*KV*NHCODE  C  17 

DO  30  K* 7, 9  C  18 

30  L(K)«L(X-l)*LA  C  19 

L  A* J  V*KV*NWCQD6  C  20 

DO  40  K» 10 , 19  C  21 

40  L I K ) «L  IK  — 1 1  *L  A  C  22 

LT0T-LI19)  C  23 

IF  (LTOT.GT.NDIM)  GO  TO  50  C  24 

RETURN  C  25 

50  WRITE  (6,601  ND I H , L  TOT  C  26 

STOP  C  27 

60  FORMAT  I  1 5H0D IMENS I ONED  8Y,I6,10H,  R EQU I  RES , [ 6 , 1H. )  C  29 

END  C  30- 
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C  LEVEL  20  CHAIN  OATE  -  7301B  12/02/3* 

SUBROUTINE  CHAIN  I M ,N , I TR , AX , AY , AZ , X , Y ,2  ,X I , ET A , ZE T A , BX , BY , BZ ,G , H,  D  1 

1XB.ZBI  D  2 

DIMENSION  Hilly  Nil),  I  TR  II )  .  AXIll,  AYI1I,  AZ  I  II  ,  XBtll,  Z  B I  lit  X  D  3 

1(11.  YIll,  2(11,  XI ( 1 1  *  ETAIll,  ZETA(l),  BX 1 1 1 ,  BY  1 1 ) ,  BZ I  II »  Gil)  D  4 

2,  Hill  D  5 

CALL  B00Y1  (H.N.ITR.AX.AY.AZ.X.Y.Z.H)  D  6 

CALL  AIREL1  I M , N , I TR , AX , A Y , AZ , X , Y , Z , H I  D  7 

CALL  STRECH  IM.N.I TR,AX,AY,AZ,X,Y,Z)  0  6 

CALL  B0OY2  I  H  ,N  ,  X  ,  Y  ,Z  ,  XI  ,  E  TA  ,  ZE  T.  A  ,  BX  ,  BY  ,  BZ  ,  H )  D  9 

CALL  SETH  I  H  ,N ,  I  TR  ,  AX  ,  A  Y,  AZ ,  X  ,  Y  ,  Z,  X I  ,  ET  A ,  ZET  A ,  BX,  BY ,  BZ,  G,  H)  0  10 

CALL  A  IRE  L  2  (  M,  N  ,G  ,H,  X6 ,  Z  B)  0  11 

CALL  USE R  V  I  M  ,N  ,  I  TR  ,  A X,  A Y  ,  AZ  •  X,  Y  ,  Z , G  ,HI  0  12 

CALL  USERS  I M ,N , I TR , AX , AY , A Z , X , Y , Z.G, HI  0  13 

RETURN  D  1* 

END  D  15- 


G  LEVEL  20  B00Y1  DATE  «  7301B  12/02/3* 

SUBROUTINE  900 Y 1  I M, N , l TR , AX , AY , AZ ,X , Y , Z  , G)  E  1 

DIMENSION  Hill,  Nil),  I TR 1 1  I ,  AXIll,  AYll),  AZ I  1 1 ,  Gill  E  2 

DIMENSION  XllV.KV.l),  Y(IV,KV,1I,  ZIIV.KV.ll  E  3 

COMMON  /CMMNA/  KL , JRU8 , IRE  ST , MN SUM, NWCODE , LXY , LXZ , L900  E  * 

COMMON  /CMMNB/  R ER , T I TL E I  20 1  , F S 

COMMON  /CMMNC/  I V , J V, KV  E  5 

00  10  K-l.NWCODE  E  6 

Ml»M(KI*ITR(Kl  E  7 

NI'NIK)  E  B 

00  10  J«1,N1  E  9 

00  10  (-1.M1  E  10 

10  READ  13,501  XI  I ,J,K) , Y(  I , J,K) ,Z< I ,J,K)  E  11 

IF(  L900.EQ.0)FS»1. 

(FIL900.EQ. 1IFS*1. 

(FIL900.EQ.2IFS*0. 

IF  (L900.EQ.OI  GO  TO  20  E  12 

IFIL900.EQ.3IG0  TO  20 

READ  (5,30)  I G ( J I  ,  J*1 » MNSUM I  E  13 

GO  TO  *0  E  1* 

20  DO  30  J«l, MNSUM  E  15 

30  G  ( J I  *  0  •  E  16 

*0  RETURN  E  17 

50  FORMAT  (6E12.0I  E  19 

END  E  20- 
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G  LEVEL 


10 


20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 


20 


AIRFL1  OAT  F  =  73018  12/02/34 


SUBROUTINE  AIRFL1  (  M,  N ,  I  T  R  ,  AX  ,  A  Y  ,  AZ,  X  ,  Y  ,  Z  ,  G) 

DIMENSION  M{  1 )  ,  *  N(  1 )  ,  ITR(l),  Axil).  Aril),  AZ I  1  I ,  Gill 
DIMENSION  XIIV.KV.ll.  YlIV.KV.il.  ZlIV.KV.l) 

COMMON  /CMMNA/  M I NGO , JRUB , I  RE  ST , MNSUH , NWCOOE , LXY , LX Z , L900 
COMMON  /CMMN8/  RER , T I TLEI 20  I , F S , AM, AG AM, AMS . F 1 , F 2 
COMMON  /CMMNC /  I  V , J V , K V ,L TOT , MTO T 
REAL  RPD/. 01 745329/ 

WRITE  16,301  TITLE 
WRITE  (6,601 

WRITE  (6,40)  NWCOOE , LXY, LXZ ,L900, MINGO, IREST, MNSUM, LTOT 
WRITE  (6,70) 

WRITE  (6,50)  FS.RER.AM.AGAM 
DO  10  K=l, NWCOOE 
WRITE  (6,801 

WRITE  (6,901  K ,M( K) ,N(KI , I TRI Kl , AX(K) ,AY I K) , AZ( K) 

WRITE  (6,100) 

M  1*M  I  K )*ITR(K) 

N  1-  N  (  K  I 
DO  10  1*1, Ml 

WRITE  (6,110)  I , ( X( I , J,KI ,Y( I ,J ,K) ,Z( I , J,K) , J  =  1 ,N1) 

WRITE  (6,120) 

LTOT*  1 
M  TOT  *  1 
RER=RER*RER 
A  GAM*. 5*1 AGAM-1.  ) 

ams=am»am 
F  2* 1 . / I  l.-AMSI 
F  1  * SQR  T  (  F  2 1 
DO  20  K* 1 , NwCODE 
AX(KI*C0SIRPD»AX(K) ) 

AY(K|*CDS(RP0*AYIKI ) 

AZ(KI=COS(RPD*AZ(K) I 

rf  turn 

FORMAT  I1H1.20A4) 

TORMAT  ( 10112) 

FORMAT  | 1P10E12.4) 

FORMAT  I 1H0, 7X, 6HNWC00E , 8 X, 3HLXY ,9X , 3MLX Z , 8X , 4HL900 , 7X , 5HMING0, BX, 
15HIREST, 5X , 5HMN  SUM ,8X,4HLT0T) 

FORMAT  ( 1H0 , T  3 , 9HF  RE  E  S  TR  M  ,  T  1  5 , 3HRE  R  ,  T27 , 2HM I  ,T  39 , 1 HK  ) 

FORMAT  ( 1H0, T4.4HPART , T 1 2 . 1 HM , T l 8 , 1 HN , T 2 3 , 3H I TR , T 39 , 2HAX , T5 1 , 2HAY , 
1T6  3, 2HAZ ) 

FORMAT  I4I6.6X.3F12.2I 

FORMAT  I  1H0.T6,  M  • ,T15,'X • , T2  7, *V* ,T39,’ Z'  ,T57,'X>  ,T69, *Y» ,T81, 'l' 
1,  T99,  •  X*  .Till  ,•  Y*  ,T123,  *Z  ■) 

FORMAT  (  I  6, 6X, 1P3E 12. 4.6X.3E 12.4 ,6X .3E12.4/ I 12X.3E12.4, 6X.3E12.4.6 
IX, 3E 12.4) ) 

FORMAT  ( 1 HA ) 

END 


F  1 
F  2 
F  3 

F  5 
F  6 

F  7 
F  8 
F  9 
F  10 
F  11 
F  12 
F  13 
F  14 
F  15 
F  16 
F  17 
F  18 
F  19 
F  20 
F  21 
F  22 
F  23 
F  24 
F  25 
F  26 
F  27 
F  28 
F  29 
F  30 
F  31 
F  32 
F  34 
F  35 
F  36 
F  37 
F  38 

F  40 
F  41 
F  42 
F  43 
F  44 
F  45 
F  46 
F  47 
F  48- 
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G  LEVEL  20  STRECH  DATE  =  73018  12/02/3* 

SUBROUTINE  STRECM  {  M,  N ,  (  T R  ,  AX  ,  A  Y  ,  AZ  ,  X  , Y  ,  Z )  G  1 

DIMENSION  Mil),  Nil).  ( TR ( 1 1 i  AXCll,  AY ( 1 ) ,  A2I1)  G  2 

DIMENSION  XI IVfKVil)  ,  YllV.KV.ll,  Z I  I  V , K  V , 1 1  G  3 

COMMON  /CMMNA/  KL , JRUB , l R EST , MNSUM, NWCODE  G  * 

COMMON  /CMMNB/  RER.Tl TLE120) ,FS, AM, AGAM, AMS, FI, F2  G  5 

COMMON  /CMMNC/  IV.JV.KV  G  6 

IF  (AM.LE.O.I  RETURN  G  7 

(F  ( AM .GT . . 9999 )  GO  TO  30  G  8 

DO  20  K=1,NWC0DE  G  9 

M 1=M(K 1*1 TR(K)  G  10 

N  1  =  N(  K  I  G  11 

DO  10  1=1, Ml  G  12 

DO  10  J  =  1 , N 1  G  13 

10  X(I,J,K|=F1«X((  ,  J , K )  G  1* 

X 1=F 1* AXI K  I  G  15 

F=1./SJRT(X1*X1*AV(K)«*2*AZ(K)«*2I  G  16 

AX(K)=F*X1  G  17 

AY(K)=F*AY(K I  G  18 

20  AZ(KI=F*AZ(K)  G  19 

RETURN  G  20 

30  WRITE  (6,40)  AM  G  21 

STOP  G  22 

40  FORMAT  ( 5M0M (  = , 1PE 12 .4 , 1 8H  GREATER  THAN  ONE) 

END  G  25- 
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G  LEVEL  20 


800Y2  DATE  ■  73018  12/02/34 


SUBROUTINE  300Y2  ( M,N ,X , V , 2 , X 1 , ET A, 2ET A , BX ,BY , BZ, G) 

DIMENSION  Mil),  Nil),  Xllllf  ETAIlli  ZETA(l),  BX(l),  BY (  1 ) ,  Bill), 
IGIII 

DIMENSION  XlIVfKVtllt  Y(IV,KV,1>,  Z(IV,KV,1) 

DIMENSION  XXI2I,  YY(2»,  22(21 

COMMON  /CMMNA/  KL  , JRUB, I  RE  ST , MN SUM, NN CODE 

COMMON  /CMMNB/  RER , T1 TLE( 20) 

COMMON  /CMMNC/  IV,JV,KV 

LAG-0 

L-0 

00  60  K- 1 1 NMCOOE 
Ml-MIKI-1 
Nl-NIKI-1 
00  60  J-l.Nl 
Jl-J*l 

DO  60  I  * l » Ml 
11-1*1 
L-L  *1 

X  I ( L !»• 25* (XI  I , J,K)*XI I , J1 ,K)+X(I 1,J1 ,K» *X  1 1  1 , J,K>> 

ETAIL l«.25*( Y(  ( , J, K)*Y( 1 , J1 ,K)*Y( I  1 , J1 ,K)*Y( I  1, J,K)  ) 

2ETAIL  >-.25*121  I , J , K) *2 (  I , J 1 , K) +2 ( I  1 , J 1 , K> *2 (  1 1 . J , K >  ) 

XX( 1 >  —  X ( I ,J1,K)-X( I ,  J  ,K) 

YY( 1 >  * Y( I , Jl.Kl-YI I , J,K> 

22( 1 >*2( I ,JltK|-2( I ,J ,  K) 

XX(2)-X( I  1 ,J,K)-X( I , J ,K) 

YY(2I»Y( IltJfKl— Y( I  ,  J  ,  K ) 

22(2)«2(  I1,J,K)-Z<  I.J.Kl 
R 1— XX( 1)**2*YY( 1)**2*ZZ( 11**2 
R2-XXI 2)**2*YY( 21**2*22(21**2 
IF  (Rl.GT.O.)  GO  TO  10 
IF  (R2.E0.0.1  GO  TO  40 
KK-1 

GO  TO  20 

10  IF  (R2.GT.0.I  GO  TO  30 

KK-2 

20  XX(KKI«X( I1,J1,K)-X(I ,J,K> 

YY(KK>-Y(  I  1, J1 ,K>-Y( I  , JfK> 

22(KK>-2< I 1.J1 *  K» -2 1 (  ,J,K> 

F»XX( KK)**2+YYI KK)**2*ZZ( KK) **2 
IF  (F. Eg. 0.1  GO  TO  40 
30  X3»YYI1>*Z2I2>-Z2U)*YYI2> 

Y3-22I 1I*XX(2I-XX( 11*22(21 
Z3-XXI 1>*YY(2>-YY( 1>*XX(2> 

F=X3*X3*Y3*Y3*23*23 
IF  (F.GT.O.>  GO  TO  50 

40  LAG-1 

WRITE  16,80)  I , J , K 

GO  TO  60 

50  F-l./SQRTIFI 

B  X(L  I  »F*X3 
BY(L I-F*Y3 
BZ( L  l»F*Z 3 
60  CONTINUE 

IF  (LAG.EQ.l)  STOP 
WRITE  (6,90)  TITLF 
WRITE  (6,100) 


H  1 

H  3 
H  4 
H  5 
H  6 
H  7 
H  8 
H  9 
H  10 
H  11 
H  12 
H  13 
H  14 
H  15 
H  16 
H  17 
H  18 
H  19 
H  20 
H  21 
H  22 
H  23 
H  24 
H  25 

H  26 
H  27 
H  28 
H  29 
H  30 
H  31 
H  32 
H  33 
H  34 
H  35 
H  36 
H  37 
H  38 
H  39 
H  40 
H  41 
H  42 
H  43 
H  44 
H  45 
H  46 
H  47 
H  48 
H  49 
H  50 
H  51 
H  52 
H  53 
H  54 
H  55 
H  56 
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G  LEVEL  20 


DATE  »  73018 


12/02/34 


L-0 

00  70  K»  1 » NMCODE 
M1-M(K)-1 
N1-N(K 1-1 
DO  70  J-l.Nl 
DO  70  I-l.Ni 
L»L*  1 

WRITE  (6,110)  I,J,K,L,XI(L),ETAIL),ZETAIL),BXIL),BY(L),BZ(L)iG(L) 
WRITE  (6,120) 

RETURN 

FORMAT  C 12HOZERO  NORMAL, 315) 

FORMAT  ( 1H1.20A4) 

FORMAT  ( 1H0,T6,1HI , T1 2 , 1H J , T1 8, 1HK , T24 , 1 HL , T3 3, 2HX I , T45 , 3HET A, T 57, 
14HZETA,T75«2HBX«T87,2HBY,T99,2HBZ»T117,7HV  DOT  N) 

FORMAT  (4I6,6X,1P3E12.4,6X,3E12.4,6X,E12.4) 

FORMAT  ( 1HA ) 

END 


AEDC-TR  73-9 


G  LEVEL  20  SETH  OATE  -  73018  12/02/34 

SUBROUTINE  SETH  (  M  ,  N,  I  TR ,  AX ,  A  Y ,  AZ  ,X ,  V  ,  Z ,  X I ,  ET  A,  2ET  A,  BX,  BY  ,  B2 ,  A  I ,  H)  I  1 
DIMENSION  Mill,  Nil),  ITR(l),  AX(l),  AY  (  1 )  ,  AZ(1),  XI(l),  E  T  A I  1 ) «  I  2 
lZETA(l),  B  X  (  1 )  f  B  Y  (  1 )  f  61111.  Aid)  I  3 

DIMENSION  HIMNSUM),  XIIV,RV,1),  YdV.KV.l),  ZlIV.KV.l)  I  4 

COMMDN  /CMMNA/  KL ,JRUB, IREST, HNSUM, NWCODE ,LXY,LXZ  I  5 

COMMON  /CHHNB/  RE k , T I TLE ( 20 ) . E S  I  6 

COMMON  /CMMNC/  I  V , J V , KV , KXY , KXZ  I  7 

IF  (IREST.GT.2I  GO  TO  170  I  8 

IF  IIREST.EO.O)  GO  TO  20  19 

DO  10  I  STAR T  =  1 , HNSUM  I  10 

READ  ( 21 . END-2 )  H,DX  I  11 

10  WRITE  (11)  H  t OX  I  12 

RfWINO  21  I  13 

GO  TO  150  I  14 

2  REWIND  21 

WRITE  (6,200)  ISTART  I  16 

GO  TO  30  I  17 

20  I  STAR  T  =  1  I  18 

30  DO  40  I " l » HNSUM  I  19 

40  A  |(  |  )«H<  I  )-FS*BX<  I  )  I  2D 

DO  140  1*1  START , HNSUM  I  21 

DO  50  J«l, HNSUM  I  22 

50  H(  J )-  0  •  I  23 

JXY*1  I  24 

JXZd  I  25 

GO  TO  80  I  26 

60  CALL  YEMY  < JXZ,M,N,ITR,AX,AY,AZ,X,Y,Z)  I  27 

GO  TO  8D  I  28 

7D  CALL  ZEMZ  <  JX  Y  ,M  ,N  ,  I  TR  ,  AX  ,  AY ,  AZ  ,X  ,Y  ,  Z  )  I  29 

80  PM=KXY*KXZ  I  30 

I T*0  I  31 

DO  IDO  Ja 1 , NWCODE  I  32 

M 1  =  H ( J  ) -  I  I  33 

M2»H(JI»ITR( J)-l  I  34 

N1»N(J)-1  I  35 

DO  100  K 1=  1  »N1  I  36 

K2»K1*1  I  37 

DO  100  L 1* 1 , H 1  I  38 

I  T*  I  T  ♦  1  I  39 

CALL  SEGH  (PM.XILl ,K1 ,J) , Y(L1 ,K1 , J) ,Z  IL1 ,K1 , J) ,XIL1,K2, J ),Y(L 1,K2,  I  4D 
1J  ),Z(L  1 ,  X  2 ,  J  )  ,XI  (  I  )  ,ETA(I  )  ,  ZE  T  A  (  I  )  ,0X  »0Y  ,DZ)  I  41 

00  90  L 2* L 1 ,H2  I  42 

L  3*L  2*1  I  43 

CALL  SEGM  (PH,X(L2.Xl.J)fY(L2.Kl,J)lZ(L2,Xl.J),X(L3fKl,J),Y(L3,Kll  I  44 
1J ) ,Z( L3.A1 . J) .XI ( I ) ,ETA ( I ) , ZETA( I ) ,XL »YL ,ZL )  I  45 

CALL  SEGM  (PMfx(L2,K2tJ),Y(L2»K2,J),Z(L2,K2,J)fX(L3,R2,J),Y(L3,K2f  I  46 
1J ) •Z(L3,K2. J) ,XI < I ) ,E TA( I ) , ZETAI I ) ,XR,YR,ZR )  I  47 

OX»DX-XL*Xk  I  48 

DY*DY-YL*VR  I  49 

90  OZ*DZ-ZL*ZR  I  50 

CALL  RAY  (PH,X(L3,R1. J) ,Y(L3,M, J) ,Z(L3,K1, J) ,AX( J) ,AY( J ) ,AZ< J) ,XI  I  51 
1<  I  ) ,ETA(  I ) , ZETAI I ) ,XL ,YL,ZL)  I  52 

CALL  RAY  (PM,X( L3.K2 , J) , Y (L3.K2 , J) , Z(L3 ,K2 ,J)  , Ax  I J) , AY  I J ) , AZ< J ) , XI  I  53 
1(  I ) ,ETA< I  ) ,ZETA ( I  )  ,XR  ,YR, ZR)  I  54 

DX*OX-XL*XR  I  55 

D  Y»OY- YL ♦ YR  I  56 
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C  LEVEL  20  SETH  DATE  -  7301B  12/02/34 

0Z»0Z-ZL*ZR  I  57 

ioo  H<iT)»Him+Dx*sxin+or*sv(i)+DZ*Bzm  i  se 

IF  IJXr.LT.OI  GO  TO  120  I  59 

IF  1JXZ.LT.0)  GO  TO  110  I  60 

IF  1LXZ.GT.0)  GO  TO  60  1  61 

IF  ILXV.GT.O)  GO  TO  70  I  62 

GO  TO  130  I  63 

110  IF  ILXV.GT.O)  GO  TO  70  I  64 

GO  TO  130  I  65 

120  IF  (JXZ.LT.O)  GO  TO  60  I  66 

130  CONTINUE  I  67 

WRITE  (11)  H  #  A I ( I)  I  68 

IF  (JRU8.E0.1)  GO  TO  160  I  69 

140  CONTINUE  I  70 

150  l-MNSUM  I  71 

160  ENO  F ILE  11  I  72 

REWIND  11  I  73 

IF  (I.LT.MNSUM)  GO  TO  180  I  74 

170  RETURN  I  75 

180  WRITE  (6,190)  I , HNSUM  I  76 

STOP  I  77 

190  FORMAT  (  1H0, 1 5, 3H  0F,I5,8H  IN  SETM)  I  79 

2 CO  FORMAT  ( 24M0 START  I NG  IN  SETM  ON  ROW, 15)  I  80 

ENO  I  81- 
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AEDC-TR-73-9 


G  LEVEL  20 


SEGM  DATE  «  73018  12/02/34 


SUBROUTINE  SEGM  ( G , AX , AY , AZ ,BX, BY ,01 ,CX , CV ,CZ ,DX . DY , D2) 

COMMON  /CMMN3/  PS 

EX*BX-AX 

FY-BY-AY 

FZ-BZ-AZ 

ES»EX*EX+EY*EY*EZ*EZ 
IF  ( ES.EQ.O.)  GO  TO  10 
FX=CX-AX 
F  V=C  V- A V 
FZ=CZ-AZ 

F  SXF  X*F  X-fF  Y*FY«-FZ*FZ 

IF  (FS.EQ.D.)  GO  TO  10 

DX*EY*FZ-EZ*F Y 

DY«EZ*FX-EX*FZ 

0Z=EX*FY-E Y*F  X 

DS«DX*DXtOY*OY*DZ*DZ 

IF  (OS.LE .ES*PS)  GO  TO  10 

HS=ICX-BX)**2*( CY-BY) **2*ICZ-BZ I **2 

F-SQRTIFSI 

H*  SQR  T (HS ) 

CT1«1 FS-HS*ES)/F 
CT2»( F  S-HS-ES) /H 
C=G*(CTl-CT2)/l 25. 1 32 74*0 S) 

OX«C*DX 
f)Y=C*OY 
DZ*C*DZ 
RETURN 
10  DX=D. 

D  Y=0. 

OZ  =  D. 

RETURN 

END 


J  1 
J  2 
J  3 
J  4 
J  5 
J  6 
J  7 
J  8 
J  9 
J  10 
J  11 
J  12 
J  13 
J  14 
J  15 
J  16 
J  17 
J  18 
J  19 
J  20 
J  21 
J  22 
J  23 
J  24 
J  25 
J  26 
J  27 
J  28 
J  29 
J  30 
J  31 
J  32- 


G  LEVEL  20 


RAY 


DATE  «  73018  12/02/34 


SUBROUTINE  RAY  I G , AX, AY , AZ , EX ,E Y , EZ , CX ,C Y ,CZ , DX , DY, 02 ) 
COMMON  /CMMNB/  PS 
F  X*  C  X—  A  X 

fy=cy-ay 

FZ=CZ-AZ 

FS»FX*FX*F Y*FY*F2*F2 

IF  IFS.EQ.O.)  GO  TO  10 

OX=EY*FZ-EZ*FY 

OY»EZ*FX-EX*FZ 

DZ=EX*FY-E Y*FX 

OS=OX*OXTOY*OY*DZ*OZ 

IF  (OS.LE.PS)  GO  TO  10 

CTl«t  EX*FX*EY*F Y*EZ*FZ) /SQRTIFS) 

OG*ICTl+l.l/(  1 2. 5663  7*0S ) 

DX»C*DX 
OY*C*DY 
OZ*C*OZ 
RE  TURN 
10  DX«0. 

0Y=0. 

OZ«D. 

P  F  TURN 
END 


K  1 
K  2 
K  3 
K  4 
K  5 
K  6 
K  7 
K  8 
K  9 
K  ID 
K  11 
K  12 
K  13 
K  14 
K.  15 
K  16 
K  17 
K  18 
K  19 
K  20 
K  21 
K  22 
K  23- 
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AEDC-TR-73-9 


G  LEVEL  20  VEMY  DATE  »  73018  12/02/34 

SUBROUTINE  YEMY  t L ,M, N, I TR , AX , AY , AZ ,X ,Y , Z )  L  1 

DIMENSION  Mill,  Nil),  1TRI1),  AXI1),  AVID,  AZI1)  L  2 

DIMENSION  XI  (  V,  KV « 1  i  «  V(jV,KV,l),  ZKV.KV.l)  L  3 

COMMON  /CMMNA/  I DUM1 4 ) , NWCODE  L  4 

COMMON  /CMMNC /  1 V , J V , KV ,KXY , KXZ  L  5 

DO  10  K«1 .NWCODE  L  6 

AYIK)  —  AYIK)  L  7 

Ml-MIK )♦( TRIK)  L  8 

Nl-NIK)  L  9 

DO  10  J— 1 ,N  1  L  10 

DO  10  1-1, Ml  L  11 

10  YI  I,  J.K)  — YI  I  ,J,K)  L  12 

L—L  L  13 

KXZ— KXZ  L  14 

RETURN  L  15 

END  L  16- 


G  LEVEL  20  ZEMZ  DATE  -  73018  12/02/34 

SUBROUTINE  ZEMZ  I L . M, N, I TR , AX , A Y , AZ t X . V . Z )  M  1 

DIMENSION  Mil),  Nil),  ITR(l),  AX(l),  AY  I  1 ) •  AZI1)  M  2 

DIMENSION  XII  V,  KV ,  1 )  ,  Y((V,KV,1),  ZKV.KV.l)  M  3 

COMMON  /CMMNA/  I DUMI4 ) ,  NWCOOE  M  4 

COMMON  /CMMNC/  1 V • J V, KV ,KX Y.KXZ  M  5 

DO  10  K«1 .NWCODE  M  6 

AZIK)  — AZIK)  M  7 

M1*MIK)KTMK)  M  8 

Nl-NIK)  M  9 

DO  10  J-l.Nl  M  10 

00  10  I-l.Ml  M  11 

10  Z I  I , J  « K )-— Z I  I « J . K)  M  12 

L—L  H  13 

KXY— KXY  M  14 

RETURN  M  15 

END  M  16- 
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AEDC-TR-73-9 


C  LEVEL  20 


A ( RE  L2 


DATE  «  73016  12/02/34 


SURROUT I NE  AIRFL2  ( M, N,G , H • A , 6 1 
DIMENSION  Mill.  Nil).  Bill 

DIMENSION  A I MNSUMI  •  G(JV,KV,1),  H(JV,KV,1) 
COMMON  /CMMNA/  ( DE T , JRUB , I  REST, MNSUM , NWCODE 
COMMON  /CMMNB/  RER.TI TLEI20) 

COMMON  /CMMNC /  IV,JV,KV 
IE  I IREST.CT.2I  GO  TO  10 
CALL  LSTSEQ  (G.H.A.BI 
(E  ( IOET.EO.O)  GO  TO  110 

10  READ  (121  A 
REMIND  12 
L*0 

DO  20  K *  1 # NMCUDE 
M1«M(KI-1 
N1*NIK 1-1 
DO  20  J*1,N1 
DO  20  ( *  1 , Ml 
L*L  +  1 

20  GII,J,KI*A(L) 

WRITE  16,1201  TITLE 
WRITE  (6,1301 
DO  30  K*l, NWCODE 
WRITE  16,1401  K 
Ml-MlKl-1 
N  1*NI  K  I  —  1 
DO  30  1*1, Ml 

30  WRITE  16, ISO)  I ,(GII,J,K>,J*1,N1) 

DO  90  K»l, NWCODE 
M1*M(KI  —  1 
N1=N(K> 

DO  90  J-l.Nl 

00  90  I  *  1 , Ml 

IF  ( I .GT. II  GO  TO  40 

Gl-O. 

GO  TO  50 

40  Gl-HI I-1,J,K| 

50  IF  I J.GT. II  GO  TO  60 

G2»0. 

GO  TO  70 

60  G2*Gl  I.J-l.KI 

70  IF  IJ.LT.N1)  GO  TO  80 

G3-0. 

GO  TO  90 

80  G3— G(I,J,K) 

90  H( I,J,KI«G1*G24G3 

WRITE  16,1701 
00  100  K*l, NWCODE 
WRITE  16,1401  K 
M  1«M(  K  1-1 
Nl'NIKI 
00  100  1*1, Ml 

100  WRITE  (6,150)  ( ,IHI I , J.KI , J*1 ,N1) 

WRITE  16,1801 
RETURN 

110  WRITE  (6,1601 

STOP 


N  1 
N  2 
N  3 
N  4 
N  5 
N  6 
N  7 
N  8 
N  9 
N  10 
N  11 
N  12 
N  13 
N  14 
N  15 
N  16 
N  17 
N  18 
N  19 
N  20 
N  21 
N  22 
N  23 
N  24 
N  25 
N  26 
N  27 
N  28 
N  29 
N  30 
N  31 
N  32 
N  33 
N  34 
N  35 
N  36 
N  37 
N  38 
N  39 
N  40 
N  41 
N  42 
N  43 
N  44 
N  45 
N  46 
N  47 
N  48 
N  49 
N  50 
N  51 
N  52 
N  53 
N  54 
N  55 
N  56 


120 

FORMAT 

I 1H1.20A4) 

130 

FORMAT 

( 1H0.T6, M • , T1 5, •GAMMAS*) 

140 

FORMAT 

( 1H0,T15, 4HP AR  T , 1 3 1 

150 

FORMAT 

( I  6, 6X.1P10E 12.4/1 12X10E12.4) I 

160 

FORMAT 

I 16H0SINGULAR  MATR1XI 

170 

FORMAT 

I 1H0.T6, M ' ,T15,»ETAS»I 

180 

FORMAT 

ENO 

I  1HAI 

N  58 
N  59 
N  60 
N  61 
N  62 
N  63 
N  64 
N  65- 
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AEDC-TR-73-9 


G  LEVEL  20 


LSVSEO 


OATE  *  73til8 


12/02/34 


SUBROUTINE  LSYSEQ  (T,A,B,L) 
DIMENSION  Till,  All),  6(1,1),  L  ( 1  > 
COMMON  /CMMNA/  I OET , JRUB, IREST, NR 
INTEGER  NS/1/, NO/10/ 

NC-NR*NS 

CALL  OAFILE  (  10, NR  ,  NC  ,  'U*  ,KL  ) 

IF  ( IREST.LT. 2)  GO  TO  10 
REAO  (22)  JLOW, (L(K) ,K*1, JLOW) 
WRITE  (6,170)  JLOW 
KK*22 
GO  TO  20 
10  JLOH-1 
K  K*  1 1 

20  00  30  J-l.NR 

REAO  (KK)  (T(K> ,K»1,NC) 

30  WRITE  (NO  *  J)  (T(K) ,  K*  1 ,  NC  ) 

PEWINO  KK 
NST-NRU-JLOW 

REAO  (NO’NST)  ( T ( K ) , K»1 , NC ) 

00  90  J- JLOW, NR 
IF  (JRUB.NE.l)  GO  TO  50 
WRITE  (12)  J ,(L(K) ,K«1,J) 

00  40  KK* 1 , NR 

REAO  ( NO  a  KK )  ( T ( K )  , K- 1 , NC ) 

40  WRITE  (12)  (T(K),K*1,NC> 

ENO  FILE  12 
PEWINO  12 
WRITE  (6,160)  J , NR 
STOP 

50  CONTINUE 

F  (NO  INO'KL-l) 

SUP-0.0 
L (  J  )- 0 

00  60  K* 1 , NR 
TLT-ABSIT(K)) 

IF  (SUP.GE.TLT)  GO  TO  60 
SUP-TLT 
L ( J )*K 

60  CONTINUE 

IF  (L(J).EQ.O)  GO  TO  150 
L  T-L  (  J  ) 

DO  70  K* 1 , NC 
70  A(K)*T(K)/T(LT) 

A( LT )»  1.0 

WRITE  ( NO ’KL  )  (  A( K ) , K-l , NC ) 

IF  (J.EQ.NR)  GO  TO  100 
FINO  (  NO  *  1 ) 

JT-NR-J 
00  90  JR-l.JT 

READ  (NO'KL)  ( T ( K ) , K* 1 , NC ) 

FINO  (NO'KL-1) 

TLT-T(LT) 

00  BO  K* 1 , NC 
80  T(K)«T(K)-TLT*A(K> 

T(LT)*Q.O 

90  WRITE  (NO'KL)  ( T( K ) ,K*1 , NC ) 


0  1 

Sil)  0  2 

0  3 

0  4 

0  5 

0  6 
0  7 

0  8 
0  9 

0  10 
0  11 
0  12 
0  13 

0  14 

0  15 

0  16 
0  17 

0  18 
0  19 

0  20 
0  21 
0  22 
0  23 

0  24 

0  25 

0  26 
0  27 

0  28 
0  29 

0  30 

0  31 

0  32 

0  33 

0  34 

0  35 

0  36 

0  37 

0  38 

0  39 

0  40 

0  41 

0  42 

0  43 

0  44 

0  45 

0  46 

0  47 

0  48 

0  49 

0  50 

0  51 

0  52 
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t  AEDC-TR-73-9 


G  LEVEL 

20  LSYSEQ 

DATE  -  73018 

12/02/34 

100 

00  110  K«1,NS 

0 

53 

KS-NR+K 

0 

54 

LS-UNR) 

0 

55 

110 

B(K,LS)>A(KS) 

0 

56 

DO  130  J-2.NR 

0 

57 

Jl-J-1 

0 

58 

JJ-NR-J+1 

0 

59 

LS-Ll JJ) 

0 

60 

READ  (ND'KL)  I At K 1 t K» 1 , NC » 

0 

61 

00  130  K»1,NS 

0 

62 

KS»NR*K 

0 

63 

S(K)«AIKS) 

0 

64 

DO  120  KK>1.J1 

0 

65 

KKK«NR-KK+1 

0 

66 

L  T»L(  KKK  1 

0 

67 

120 

SIK)*S(K)-A<LT)*B(K,LT) 

0 

68 

130 

8IK,LSI»S(K) 

0 

69 

FIND  IN0*11 

DO  140  K= 1 »NS 

0 

70 

140 

WRITE  (12)  ( BIKtKK) »KK«1,NR) 

0 

71 

END  F ILE  12 

0 

72 

REWIND  12 

0 

73 

IDET-1 

0 

74 

RETURN 

0 

75 

150 

I DET*  0 

0 

76 

RETURN 

0 

77 

160 

FORMAT  (1H0,I5,3H  0F,I5,10H  IN 

LSYSEQ) 

0 

79 

170 

FORMAT  I 26H0START I NG  IN  LSYSEQ 

ON  ROW  1 1 5 ) 

0 

80 

END 

0 

81- 

113 


AEDC-TR-73-9 


G  LEVEL  20  USERV  DATE  »  73018  12/02/34 

SUBROUTINE  USERV  (  M , N ,  I  TR  , AX ,  AY  ,  AZ ,X  , V  ,  2 , G,  HI  P  1 

DIMENSION  Mil),  Nil*.  ITR(l),  AX(l),  AYI1)  ,  AZ(l),  XU).  Y(L),  2(1  P  2 

1).  GUI,  HU)  P  3 

COMMON  /CMMN8/  RER . TI TLEI 201  P  4 

PEAL  SGN  /-l./ 

IU-13  P  5 

READ  (5,50)  NV  P  6 

IF  (NV.EQ.O)  GO  TO  30  P  7 

IPAGE-0  P  8 

ILINE«62  P  9 

00  20  J»1,NV  P  10 

READ  15.60)  XV.YV.2V  P  11 

CALL  VELOCY  1  M,  N,  I  TR  ,  AX  ,  AY ,  A2  ,X ,  Y  ,2 ,  G  ,H,X  V ,  Y  V ,  2  V.U,  V ,  M  )  P  12 

IF  ( ILINE.LT.62)  GO  TO  10  P  13 

1PAGE«IPAGE*1  P  14 

ILINF«9  P  15 

WRITE  (6,70)  TITLE, IPAGE  P  16 

CALL  VELLA8  P  17 

10  ILINF«ILINE*1  P  18 

20  CALL  VELOUT  ( X V, Y V ,2 V.U ,V , W, I U)  P  19 

WRITE  (6,80)  P  20 

30  IF  (IU.EO.O)  GO  TO  40  P  21 

WRITE  (IU)  (SGN,J»1,12)  P  22 

END  FILE  IU  P  23 

PEWIND  IU  P  24 

40  RETURN  P  25 

50  FORMAT  1 241 3)  P  27 

60  FORMAT  (6E12.0)  P  28 

70  FORMAT  ( 1H1 , 20A4 , 5X , • VELOC I T I ES * , 5X , » PAGE* , 1 3 )  P  29 

80  FORMAT  (1HA)  P  30 

END  P  31- 


G  LEVEL  20  USERS  DATE  -  73018  12/02/34 

SUBROUTINE  USERS  ( M ,N , I TR , AX , AY , A2,X , Y , 2,G,  H)  Q  1 

DIMENSION  Mil),  Nlll,  I  TR  ( 1  I  ,  AXU),  AY  (1 )  ,  AZ(  1 )  ,  Xlll,  Y(l),  211  0  2 

1),  G(  1),  HU)  0  3 

REAL  SGN  /-!./ 

I U* 14  04 

READ  (5,40)  NS  05 

IF  (NS.EO.O)  GO  TO  20  06 

DO  10  J«1,NS  Q  7 

PEAD  (5,50)  XO, YO.ZO.OSO, AN , AM, NX,XN, XX , YN , YX , 2N, 2X  0  8 

CALL  STREAM  (  M , N ,  I  TR ,  AX  ,  AY , AZ  ,X  , Y  ,Z,G ,H,X0 ,  YO  , 20 ,  DSO ,  AN  ,  AM,  NX  ,XN,  X  Q  9 

1X.YN, YX.ZN.ZX, I U)  Q  10 

IF  (IU.EO.O)  GO  TO  10  0  11 

WRITE  (IU)  (SGN,K*1,12)  0  12 

10  CONTINUE  Q  13 

20  IF  (IU.EO.O)  GO  TO  30  0  14 

END  FILE  IU  0  15 

REWIND  IU  0  16 

30  RETURN  Q  17 

40  FORMAT  (2413)  Q  19 

50  FORMAT  ( 6E 12. 0, I 8/6E12. 0)  Q  20 

END  Q  21- 
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AEDC  TR  73-9 


G  LEVEL  20  VELOCY  OATE  *  73018  12/02/34 

SUBROUTINE  VELOCY  ( M , N, 1 T R , AX .A Y , AZ *X ,Y , Z ,G .H , A B , YO , 20 , T U , T V , T W )  R  1 

0 1  ME  NS  ION  Mil),  Nil),  I TR ( 1 ) ,  AXllI,  AY (11 .  AZtl>  R  2 

(DIMENSION  X(IV.KV.l).  Y(IV,KV,1I*  Z(IV,KV,1>,  G(JV,KV,1I,  HIJV.KV,  R  3 

11  )  R  4 

COMMON  /CMMNA/  KL , JRUB, I R EST , MNSUM.NWCOOE , LXY , LXZ  R  5 

COMMON  /CMMN B/  RER , T 1 TLE ( 20 1 ,F S , AM, AG AM, AMS » F 1 , F 2  R  6 

COMMON  /C  MMNC /  I  V ,JV ,KV ,KXY ,KXZ  R  7 

XO=F 1*  XB  R  B 

TU=FS  R  9 

TV®0.  R  10 

TN=0.  R  11 

JXY=1  R  12 

JXZ=1  R  13 

GO  TO  30  R  14 

10  CALL  YEMY  ( JXZ , M , N , I TR  ,  AX , A Y  ,  AZ  ,X  ,Y , Z I  R  15 

GO  TO  30  R  16 

20  CALL  ZEMZ  ( JXY , M, N , I T R , AX , A Y , AZ , X , Y , Z I  R  17 

30  JPM-KXY*KXZ  R  IB 

00  80  L  = 1 .NWCOOE  R  19 

N  1*N(  LI  —  1  R  20 

N  2=N  (  L  )  R  21 

M1=M(L)-1  R  22 

M2«M(  L  1*1  TR(L  1-1  R  23 

00  40  Kl»l,Nl  R  24 

K  2-X 1  +  1  R  25 

00  40  J1=1,M1  R  26 

GAM=JPM*G( J1 ,K  1 ,L)  R  27 

CALL  SEGM  ( GAM , XI J 1 , K 1 , L ) , Y ( J 1 , K 1 ,L1 , Z ( J 1 , X 1 , L I , X ( J 1 , K2 , L ) , Y I J 1 , K 2  R  28 

1,L),Z(J1,K2,L) ,XO,YO,ZO,U,V,W)  R  29 

T  U®  TU  +  U  R  30 

T  V®  T  V  + V  R  31 

40  Tw=TW+W  R  32 

00  70  K 1* 1 ,N2  R  33 

00  60  J 1® 1 ,M2  R  34 

J  2= J 1+ 1  R  35 

IF  (J1.GT.M1I  GO  TO  50  R  36 

GAM®JPM*H(J 1,K1  ,L)  R  37 

50  CALL  SEGM  ( GAM , X ( J 1 , K 1 , L )  . Y ( J 1 , K 1 , LI  , Z IJ 1 , K 1 , L 1 , X  I J2 , K1 , L I , Y ( J2 , K  1  R  38 

l.Ll.ZI J2,K1,LI,X0,Y0,Z0,U,V,W)  R  39 

TU®TU+U  R  40 

TV®TV+V  R  41 

60  TH=TW+W  R  42 

CALL  RAY  (GAM, XI J2.K1 ,L I ,YC J2 ,K1,LI ,ZCJ2 ,K1,U ,AX(L) , AV (LI, AZ(L 1,X  R  43 

10, YO.ZO.U, V,WI  R  44 

T  U=  TU  +  U  R  45 

T  V®  T  V+  V  R  46 

70  T  W®  T  W  +  W  R  47 

80  CONTINUE  R  48 

IF  (JXY.LT.OI  GO  TO  100  R  49 

IF  ( JXZ.LT.OI  GO  TO  90  R  50 

IF  (LXZ.GT.OI  GO  TO  10  R  51 

IF  (LXY.GT.OI  GO  TO  20  R  52 

GO  TO  110  R  53 

90  IF  (LXY.GT.OI  GO  TO  20  R  54 

GO  TO  110  R  55 

100  IF  (JXZ.LT.OI  GO  TO  10  R  56 

110  TU«1.+F2*(  TU-1.  1  R  57 

TV«Fl*TV  R  58 

TW®Fl*TW  R  59 

RETURN  R  60 

ENO  R  61- 
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AEDC-TR-73-9 


G  LEVEL  20  VELOUT  DATE  ■  73018  12/02/34 

SUBROUTINE  VELOUT  ( XD , YO, 20 , U , V , W, I U J  S  1 

COMMON  /CMNN8/  R ER  ,  TI  TLE (  20)  , FS  ,  AM,  AG  IH,  AMS  S  2 

REAL  DPR/57.29578/ 

TUS«U*U*  V*V*W*W  S  3 

TOSQRT(TUS)  S  4 

ACH-1.*AG1H*AMS*( l.-TUS)  S  5 

IF  (ACH.LE.O.)  GO  TD  10  S  6 

ACH»AM*TU/SQRT( ACH)  S  7 

10  CONTINUE  S  8 

T V*OPR*AT  AN2 ( V, U)  S  9 

TW*DPR*ATAN2( W,U)  S  10 

PC»1.-TUS  S  11 

D* ACH-AM  S  12 

IF  ( IU.EO.O)  GD  TO  20  S  13 

WRITE  (IUI  XD,YD,Z0,U,V,W,TU,ACH,TV,TW,PC,0  S  14 

20  WRITE  (6,401  XD  ,  YD ,  2D  ,U ,  V ,  W,  TU,  ACH ,  T  V  ,T  W  ,  PC  ,  D  S  15 

RETURN  S  16 

ENTRY  VELLAB  S  17 

WR(TE  (6,30)  S  18 

RETURN  S  19 

30  FORMAT  ( 1H0.T4, •X»,T16,»Y,,T28, •Z,,T41,*U»,T51,,V»,T61, »W*,T70, • |V  S  21 

II  *,T81,  'MST88,  *ATN(  V/U)  *  ,T97  ,*  AT  N(  W/U)  •,  T  1 08  ,*  CP  * ,  T  1 20 ,  *M-M  I »  )  S  22 

40  FORMAT  ( 1P3E12.4,0P5F10.5,F8.2,F9.2,1PE12.3,E11.3)  S  23 

END  S  24- 
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G  LEVFL  20  STREAM  DATE  =  73018  12/02/ 34 

S  U8RUUT I NE  STREAM  (  Ml  ,M2  .  M3  .  A 1 ,  A2  t  A3  ,  A4  ,  A5  ,  A6  ,  A7  ,  A8  ,  X  0,  V  0.  Z  0  ,  OS  0,  A  T  1 
IN, AX, NX, XN, XX,YN,YX,ZN,ZX , X  U )  T  2 

OIMENSION  Mill),  M2  ( 1  )  ,  M3(l),  Alll)  ,  A2(l),  A3(l),  Adi),  A5(l),  T  3 
1A611  )  .  A7(  1)  ,  A8(  1  )  T  A 

COMMON  /CMMN8/  RER , TI TLE ( 20 )  T  5 

PEAL  RPD/. 01745329/ 

C  DC  TOOD  8/11/71  T  6 

REAL  CGOEI 8) /*N>NX* , ' XCXN1 , 'X>XX* , »Y<YN' , *Y>YX* , • Z<ZN* , * Z>ZX» , • 0<D  T  7 


IN*/ 

L=  1 

IPAGE=0 
1L  [NE  =  62 
D SX=  A8  S(  0 SO ) 

DSN=.001*0SX 

CX=COS(RPD*AN) 

CN=COS(RPD*AX) 

D  S*D  SO 
N  =  0 
X=XO 
y=yo 
z=zo 

10  CALL  VELOCY  < Ml , M2 , M3 , A  1 , A2 , A  3 , A4 , A5 
[F  ( [LINE. LT. 62)  GO  TO  30 
IPAGE= [ P AGE* 1 
l  L  [NE  =  9 

WRITE  (6,140)  TITLE, [PAGE 
IF  (  1  PAGE .GT. 1)  GO  TO  20 
WRITE  (6,150)  DSX,AN, AX 
I L [NE - [LI NE  *2 

20  CALL  VELLAB 

30  CALL  VELOUT  ( X , Y »Z ,U» V , W, I U) 

1L[NF=IL[NE*1 
IF  (N.GE.NXI  GO  TO  130 
IF  (X.LE.XN)  GO  TO  120 
IF  (X.GE.XX)  GO  TO  110 
[F  (Y.LE.YN)  GO  TO  100 
[F  (Y.GE.YX)  GO  TO  90 
[F  (Z.LE.ZN)  GO  TO  80 
IF  (Z.GE.ZX)  GO  TO  70 
T=SQRT(U*U*V*V*W*W) 

40  F=OS/T 

X1=X*F*U 
Y1=Y*F*V 
Z  1  =  Z*F*W 

CALL  VELOCY  (Ml ,M2,M3,A1,A2,A3,A4,A5 
T 1=SQRT( U1*01*V1*V1*W1*W1 ) 
C=(U*U1*V*V1*W*W1)/(T*T1) 

IF  (C.GE.CN)  GO  TO  50 
0S= . 75*0  S 

IF  (ASS(OS).LE.DSN)  GO  TO  60 
GO  TO  40 

50  F=.5*DS 

X=X*F*( U/T*U1/T1) 

Y* Y*F* ( V/T+V1/T1) 

Z=Z*F* ( W/T+Wl/T 1 ) 

N*N*  1 


T  8 
T  9 
T  10 
T  11 
T  12 
T  13 
T  14 
T  15 
T  16 
T  17 
T  18 
T  19 
T  20 

A6,A7,A8,X,Y,Z,U,V,W)  T  21 

T  22 
T  23 
T  24 
T  25 
T  26 
T  27 
T  28 
T  29 
T  30 
T  31 
T  32 
T  33 
T  34 
T  35 
T  36 
T  37 
T  38 
T  39 
T  40 
T  41 
T  42 
T  43 

A6,A7,A8,X1,Y1,Z1,U1,V1,W1)  T  44 

T  45 
T  46 
T  47 
T  48 
T  49 
T  50 
T  51 
T  52 
T  53 
T  '  54 

T  55 
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G  LEVEL  20  STREAM  DATE  »  73018  12/02/39 


IF  (C.GT.CM  DS-SIGN(AMINl(DSX, 1.5*ABS(DS) 1 i DS ) 

T 

56 

GO  TO  10 

T 

57 

60 

l=lm 

T 

58 

70 

L  *L  +  1 

T 

59 

BO 

L  =  L  ♦  1 

T 

60 

90 

L=L*1 

T 

61 

100 

L  =  L*1 

T 

62 

110 

L  =  L  1 

T 

63 

120 

L  =  L  ♦  1 

T 

69 

130 

WRITE  (6,1601  CODEILl 

T 

65 

RF  TURN 

T 

66 

190 

FORMAT  (  1H1 ,20A9, 5 X, 'STREAMLINE • , 5X , • PAG E • ,  1 3 ) 

T 

68 

ISO 

FORMAT  <  1H0, 1PE 12.9.19H  MAX  STEP  SI Z E ,0PF20 . 3 , 2 1H  DEG  MIN  ANGLE  CH 

T 

69 

1ANGE,F20.3,21H  DFG  MAX  ANGLE  CHANGE! 

T 

70 

160 

FORMAT  ( 8H0  *****  ,A9,6H  ****♦! 

T 

71 

END 

T 

72- 

1 18 
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